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ABSTRACT 
Title of Thesis: Pointwise Bounds in Parabolic and 
Elliptic Partial Differential Equations 


Lieutenant Fred James Bellar, Jr., United States Navy, 
Doctor of Philosophy, 1961 


Thesis directed by: L. EB. Payne, Research Professor, 
Institute for Fluid Dynamics and 
Applied Mathematics. 

A method is presented for obtaining explicit upper 
and lower pointwise bounds for the solutions of rather 
general interior boundary value problems. The differential 
equations associated with these problems are of the elliptic 
type in certain sections of the thesis while both linear and 
non-linear parabolic problems are the subject of investi- 
gation in other sections. The bounds are in terms of the 
integrals of the squares of known functions and hence, in 
the linear case, improvement is possible using the Rayleigh- 
Ritz technique. 
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CHAPTER I 
INTRODUCTION 


In [er] Payne and Weinbeger presented a method for 
determining pointwise bounds for the solution of Dirichlet 
and mixed type boundary value problems for arbitrary second 
order elliptic equations and rather general domains. In 
the present paper the writer will obtain similar results in 
the parabolic case for Dirichlet and mixed type boundary 
value problems where the given data is a function of the 
domain for some problems while for others the data may depend 
on both the solution and the domain of the problem. 

The methods of this paper yield error estimates in 
terms of the integral of the square of known functions and 
thus the results lend themselves to improvement by means 
of the Rayleigh-Ritz technique when the given data is a 
function of domain (the linear case). In such cases it is 
an unfortunate fact that, even if the set of functions to 
which ons applies the Rayleigh-Ritz technique is complete 
in the function space of the solution, the actual solution 
of the problem can not be obtained using our present-day 
computers due to the machine creation of large round off 
errors. It is quite possible however that future 


generations of these computing devices will eliminate 





such errors and thus make the techniques of this paper of 
greater practical importance. In this regard we note that 
the current electronic computing machines have made of 
practical interest problem solving techniques which were 
previously only of theoretical interest since feasible uses 
were not apparent at the time of their evolution. 

For many cf the problems considered in this thesis, 
the existence of a solution has not been established under 
the general conditions necessary for the existence of the 
pointwise bound. When ee addition requirements or 
modifications are placed on the boundary value problem, by 
referring to the literature ae [5]; 9] to [15] [22 | and 
[23] we nots that the existence of the solution is guaran- 
teed. 

As was pointed out above, the estimates of this paper 
are in terms of the integrals of known functions. For 
another method of estimating the solution function of 
boundary value problems, the reader is referred to the 
excellent works of Fichera 16 | to [8 | where the solution 
function is bounded in terms of the maximum norm. Using 
the results of Fichera as a starting point, Lieberstein 19] 
derived a numerical procedure by which error bounds were 
obtained using a digital computer. 

In Chapter II we obtain some preliminary results 
which shall form the basis of many of the conclusions of 


the following chapters. The identity of Hormander [16] . 
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which was used by Payne and Weinberger [21] , is modified in 
order to allow its use in a Dirichlet type parabolic problem 
and functions analogous to (.1) and (1.18) of [21] are 
introduced and investigated. Chapter III concerns itself 
with non-linear normally parabolic problems of both the 
Dirichlet and mixed types. For the purposes of this thesis 
a normally parabolic differential operator has the form 
ieee 3(w) => 2, (ats bu.) - Zz ow 
i Q ox oy 


~ 
- 


where the matrix a1J is positive definite and 4 is strictly 
positive. We shall see that the desired estimates are 
obtained by multiplying the differential operator by a 
parametrix function and then using Green's theorem to 
generate the solution at the point in question. The unknown 
terms in the resulting expression are estimated as a con- 
sequence of the results of Chapter II and by further use of 
the divergence theorem. In Chapter IV the problems considered 
are degenerate in character, thus the matrix alJ becomes 
positive semi-definite for some problems while the function 

% need only be non-negative for others. The next to the last 
chapter of the paper outlines a method by which solution 
bounds may be generated for a class of elliptic problems 
which were not considered by Payne and Weinberger in their 
investigation of such problems 21] > In Chapter VI the 
results of the thesis are applied to obtain bounds for the 


solution of a specific boundary value problem. 





CHAPTER IT 
PRELIMINARY RESULTS 


1. General Definitions 

Let V be an open (N + 1) dimensional domain of real 
variables (x, y} = (xt, x-... xN, y) which is bounded by 
two hyperplanes, y = 0 and y = Yeu and a surface S lying 
between these planes. As usual we denote the closure of V 
by V and assume that the divergence theorem holds in V. The 
reader is referred to Chapter IV of [18] for information 
relating to the most general regions in which this theorem 


is valid. Let 


2.1 OD {x,y)) y =» (xsVIEV - 3} 
then 
2.2 V = VUD(o)UD(y,)US = VUD(o)UD(y, US, 


On the boundary of V the Euclidean outward normal is denoted 


by 
2.3 (ns 2m) = (nj » Ny + + + Ny no) 
where 
N 
2 2 
2s n- +> n>, =l, 
l 5 tani 
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Since V is a region in which the divergence theorem is valida, 
it is clear that 2.3 is (uniquely) defined almost everywhere 
on the boundary of V. In addition, unless otherwise indicated, 


we require that S be such that 


N 
2.5 min > ny = m,>o 
ot 


for all (x,y)€S at which the normal vector is defined. 

A function Z = Z(x.y)is said to be piecewise continuous 
in V if V may be divided into a finite number of subregions 
in each of which the function is continuous and has a finite 
limit as the boundary is approached. The above is the usual 
definition of piecewise continuity; the following non- 
standard definitions shall also be gmployed since they 
greatly simplify the presentation of the problem. When a 
function 2(x,y) is said to be continuously (piecewise con- 
tinuously) differentiable in x for (x,yJ€V then % is con- 
tinuous on D(y) for o<y<y, and is piecewise continuous in V 
while the 22:, if=1, 2... N, are continuous (piecewise 
continuous) on D(y) for oy<y, and piecewise continuous in 
V. Thus if 4 is twice continuously differentiable in x and 
piecewise continuously differentiable in y for (x,y)¢€V then 
% is continuous in V, the 22. i=1, 2. .e WN are con- 
tinuous on D(y) for o<y<y, and piecewise continuous in V, 
the 2% », ij =1, 2. =. . N, are continuous on D(y) 


axtyxé 
for %nKYCLY, and piecewise continuous in V, and finally <= 


OZ 
sy 








is piecewise continuous over the entire volume. 


ce. The Parametrix of the Parabolic Problem 


For (x,y)€V let the following differential operator be 
defined 
: 2 
226 TY) = (are ¥, ),5 + Oy (ZY) 


where the symbol ,; indicates partial differentiation with 


¥, the summation convention is used throughout 


respect to x 
and the components of the symmetric matrix aij = alj(x,y) 
are piecewise continuously differentiable in x for (x,y)€V 


and are such that 


N N 
2a 3 2 
Zt ofa,2 $5 Fath bea, Dae 
z i 
for any numbere (§,,.. - 8 ‘«) and all(x,y)EV with equality 


on the left holding if and only if ®,; = 0 for i=l, 2...N. 
The function 4% =%(x,y) is piecewise continuously differentiable 
in y for (x,y}¢V. For the point p = (x,,y9)€D(y,) , at which 

a bound for the solution of certain boundary value problems 
shall later be obtained, we assume the functions aid and % 
satisfy a Lipschitz condition for (x,y)€V; that is, we 


Suppose there exist positive numbers ald and A such that 


; < N r 
2.8 | ati(x,y) - atJ(p)| {y,-9)? + S (xtanky Ft att 
_ 


l 
N a 
209 | Z(x,y)- Z(p) {y,-9)? + > (xt-xb)} SA 














for (x,y)é€V. 
The parametrix function J ,, is defined with respect 


to the point p by 
i Neo2 


2.10 Bi=|hTT ty,-) |? 200)" 2 VT) ox Peorp*Qiy-n)4 


for 
Cell p= (X55 TF JED(Y,)» (x,y)éV, and Z(p)>o 


where : 


2.12 po = a,,(p) (x7 = xf) (x° ~ x8) 


and a(p) is the determinant of the inverse matrix a5 3(P)- 
On the surface Diy) “ {p} the parametrix and its deriva- 
tives are definsd as the limits of these functions as 


(x,y) 9 (x,y,) with (x,y)€V ana Xena 


Theorem 2.1. The parametrix a has the following 
properties: | 

a) At all points (x, y)EV - {p}, the parametrix is 
twice continuously differentiable. 


b) The function 
Ran 
2 =} 2| ee C 
Pals (yo°y) exp k (yee s(p)p"| T(Y.) 


is integrable over V for S and 2 satisfying 


2.1 Seo 5 arid O<A., 





c) For any bounded function Y which is continuous 


on VUD(y,) the parametrix is such that 


2.15 Z = 
a, J. TZ as = YE)= Yoo). 
Y < "bo 
Proof: From the definition of Up» it is obvious that part 
(a) is satisfied... To establish part (b) of the theorem we 
note that for c,,i = 2, 3... . appropriately defined (the 


constant c, is reserved for later use) 


2.16 J Nege a orp (Sy 4) 2p e | J F (Xp) av 
= hg a 
<5 [Seregy 2 2S Oss Hy DEED Gegeg > 
+( a2 Ys < <2cpa.? ‘a, LP tne XHeyogh +eepla’ fa (OG xg) » 


+ 66x dnt) - BEE C4(yermy> “howp|-a- Sep tys sy |dv 


ic c ve -¢- zs 
att , } 1h cya) 


4 BALE 


+04 (Yon) Ug as AAP at) = 266) +5 %i5 God) 
tH a) Was ie DGycyfy+@ (Sgata, GE) a. Pray) 


» (xx) + 2(p) 4 ha (p) (Ry eh Comm lrery = h 22 (cex*) 





“Keene (ge a gerd erp - S)24 C4 444) [du durdy 
‘to é + Ne a - 

< ih x, Sf yeast e+ cp (eas 

+on (yon) go 44A +, — 4 = se ane 


+ De4a|- _ go" By Jind co dig 


where 


im ae 
1 
2.17 re = s (x = x) 


; ; , 
2.18 gy] (x¥-x8)+8 4 (y,-y)=a0) (9) -8 J(p) iej=Hi,2 . 


with 
aa N 
i Lye 22 i 
on = (Mr! [arate Loe) {av g)-% “oh AS RF, oN 
29 
‘t = (y.=4)/ (ym) Sob) fai ait } 7 
Sm, = O07 “S l(yens) + 0 Cx, M4) — Cp) Wh Sly eres 
and 


Tayi: < = : 
2220 h, (x “xO ) ++ hoy, y) 3(x,y) B(p) 


with , ei 
h, = (x*- qo et x i {2 (x, 4) “ee (De L> { MN 
Baek 


A= sles Sob FT feqyy-ecp} 
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Now since the atJ and % satisfy 2.8 and 2.9, the functions 
gids a, h, and he are bounded over V and thus the constants 
on the right of 2.16 exist. By integration by parts with 
respect to r and then y, it is easy to establish (using 2.1h) 
that each of the integrals on the right of 2.16 is bounded 
and thus part (b) of the theorem is proved. 

To establish the validity of 2.15 we suppose that 


fo is chosen so that for O<E<E, we have 


. € 
aoe Poe i yzsees = eee) fuer als 
m 4 >to Ux, 0 pF 


“4? yo 04) 

(x) EV GEV 

ne os PRY, ds 
4? be dv)-fap|erromary — & 
(64) EV 


E 
2B Ow ° 
& 4) EV 


Since 2.22 is true for every € such that OS e<é, and since 


Z and / are continuous at peD(y,), we have that 


es) ee 2%, ds 
4 yo DCs) 
Ours) EV 


e 
= Din Lom J ee Ye ALG 
&—0O 7G 0 Wy) O 
O&, S)EV 








EL 


= Dm Dime VICE 4) J f lercy-p|e peat ee 2p A #454) Y ldnduo} 


£0 4 "Go Wo 


ran YC < MCEM)S map 
{(,¥)| 2 ‘bi zy ese “x pe e*} 


M NV 
= a a. . -* wW- a =| 
Cp) Z(p) (417) Pe ie 4 6 Ors) i, tp repr (Heyer) lard w 
S 2 ¥- bed ; p : N-2 - Y + —) <j 
vp) 2¢p)* "(4r1) * Line bie ere Caen) oep| HE Heer] ab 


teh So rad 
Pt J iw -4) 2" Lu 


+ bin a. § i 2 (w- Zn “hy oe upc 22 (4045-4) [ern doo 


= vepaca? Cry = { hewn als S Sawaya yg?" 


: ayplecpya?(4 AS “) | an o w } 
= vipecp)® Cam) ® {Bice Dom § C22 (BYEt)as 


Lepr Moto Wy © 


. (4,- my) z72 Lyp Lacey aPHCy,-re)) | dn ol he 


Hence if N is even we obtain 


2.21 


47"go  DCY) 
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é nk 
= Vp) 2p) (4M) F Lin eA } fae (3)! gor 4) 


| ae | 2(p) Near) | Cie ad uJ) 


=~ Y(p) (47) “hors bom jt (52) agoleme ears) [dw 


+ YC Yelm Bim 64°% (%S4)/ 
en E>d bone 19 CE) J4v} 


=k vip? P(#) W, 


where in the last equation it is clear that OR is the 
surface of the Nedimensional unit sphere and [’ is the 
well known gamma function. On the other hand 


for N odd we write using 2.23 


2-25 omy SS V8% pds = Hp) Bp)? (4M) ili J fi 2, 
tye OC) 0 740M 
{ >) 3 Pe ame z} (aot) Fonph ete Gerd) |e lw, 


To evaluate the integral on the right of the last equation 
it is necessary to make the following change of variables. 


Let 


\"l 
2.26 = = Zp) Vey) 
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then 
2.27 pe = (4,-4) 
and 
g ou ~| 
2.28 oe s (4-4) = spol 274-9) | wn 
2() 6° 
are — 
= Leow .2(p) *6 ape) de 
yy, F a 


i} 


z(pye fo sip (-8) Le 


0 


zcpy* P(E), 


The value of the above integral is obtained by direct 


calculation. Substitution of 2.28 into 2.25 now yields 


a 
2:29 Bioy ff 42% pds= pen *(ENE)..20Z) Wy 


_v 
= 2 ¥(p—)(m) * (EZ) CW), 


where use has been made of the well-known recurrence 
relation for the gamma function. The desired result is 


now obtained by noting from page 223 of [2 | 


2.30  () = am *r(3)] 











Wy 
and substitution of this expression into 2.2) and 2.29 yields 


: ie YZ As & Xo 
sate ‘4 e I fe #5 = Wee) 


for N either even cr odd. 


3. An Auxiliary Function for the Parabolic Problem 

Here our aim is to construct a sufficiently smooth 
function 7 went Gh will be employed in the following chapters 
to obtain solution bounds for the parabolic problem. The 


desired function Gr will be defined so that 


2.32 “Fa < 00 for "e,7)ev 
and 
2.33 J (4) <2 for (x,y) EV 


where J satisfies the conditions of the previous section. 

In addition to the above, it is necessary that <P wave a 
singularity at peD(y_)s which is of a slightly lower order 
than the singularity of the parametrix ¥ p at p € Dy.) - 

The desired function ig constructed by first considering for 


N71 

Ee seat = 
2 ei CE = Cr 4) engl taw-ecpye*(an(y,--9)) | 
and computing for (x,y)é€V 


ae : cf 
ara! IG) = (n49- ay © lh, al ND ECp)A, CE)EE") C4y-r4) 











15 
| Ly 5 es 
“EN oo ai) (46 ‘4)+ iene (@v-) z (psa LOGE) a? (xex, ) 


CB) ecysr- gen Deece® + ZB (yea) 
Fae Lan “DECH) OX BN Cpe: 4) | 

= Gor ~"3) = E Wy (2N-) 2 (j0)4, Leet Sede) - BC) (ye4) 
ena ag® oe 5 ‘a, Ce) (4 = M4) is ia(RN- d2(n4 a ig Xe a,,C)(45"8) 
(EDL orga EDK sat + 5 Gord) 

-(4n) “(2w-D2Cp) + (4) “An-)recpye at (x" Ls chika ces) 
+ (yn) “(amet Yecp) 4 i dy {) &Ex*) 2. Cp) ~~) (4.- A4) 
- (8N) (an-/) 26p) £,ce=ng Ya (BN) 2¢p) A 0° Cy.-a4)]- 


“Ho I-(aN -!) z(40)0°( 6N(4,- 6) | 


with gid, ays nh, and hb defined by 2.19 and 2.21. 

For (x,y)€V and sufficiently close to p, the second term 
in the brackets on the right of the last equation dominates 
the first and third through seventh terms while the eighth 
dominates the remaining terms. Due to the boundedness of 
the g's and hts and the piecewise continuity of a‘! in V; 
it follows that there exist positive numbers 8 4 ane Po such 


that 
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J(#) <0 (x, 4) € ix4|¢ Pc » OF 40-9 S Sor 
2 36 
IG) =O Gry) EL x, 44 0 % > OSES m4 € Sa 


Let Ss be such that 

2.37 Fyre () ged) 2 - QNDE) A <O 
for all (x,y) e€ x,y |(xsy}€V, o¢ ionse Sx + » We set 
2.38 5, = min | Se, = 


and then define the desired auxiliary function as follows in 


terms of non-negative constants % and ke 


sie _N-/ 2 “4 
“4 +2 (4y-g) F aya t-(ew-i) 2cp) A’ emarg)) $f 







OsyyEVy O<A - 4 £85 


-s (w-1+2h) 


2.39 Gd =g, (psd) Compe {-en-i) 2p)" [8m (gorg)l 
+O buf {- (2-1) 2() A. [zn (doa) %) 


G4) EV 5 80 < Yo-Ay 





Le 


The constants & and g appearing in the last equation shall 
be explicitly chosen in the proof of the theorem which 
follows this paragraph. On the surface Dlyol-t{p5; 

function 4? and its derivatives are defined as the limits 
of these functions as (x,y) approaches (Xs y,) with (x,y) éV 


and X FX,3 thus we have 


20 x,y.) = Vim BH (xy) K FX, 


Yo 
oli). (x54) = Yim A (x5y) X# Ky 
Iio 
2.2 ys = LinF, , ;(%29) x £ Xo 
J Yo 
and 
SF (%sFo) 9 OP eT) 
es ——= 
43 Oy gs roy ‘ P%o 


Theorem 2.2. The function Bg and t- appropriately 
fixed) is twice piecewise continuously differentiable in x 
and piecewise continuously differentiable in y for (x,y)€V 


except at p = (x } and is such that 


03 Vo 


2 hh hoo ,IF)< (xy) eV 


and 
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Ane Ho, WF) <o (x, yl€jxsy|(x,y) € Vox Ex y,-y} 


where the coefficients of J satisfy the requirements of 


section 2 of this chapter and, in addition, the condition 


2.6 min 2 = m5 >0 
V 


Proot: The function “> obviously satisfies the continuity 
and differentiability requirements of the theorem. 

Our aim shall now be to establish 2.ll and 2.5 for 
properiy chosen constants & and & . From 2.36, 2.37 and 
the definition of “” it is clear that the desired in- 
equalities are satisfied for all “2 o where o<yo-y ¢ 8) 
and P2P >: For [Pe and o< yo-y ¢ 8. we compute by 


means of 2.39 


Nts 


2647 TQ) = IF) +e (g-a) 7 | 5 (ty) 
+(%)yo- az - @r-Dapad(eny! 2 }- 


+ a1p2 1 (2N-) B(p)A*(BN (yp) 
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For (x,y) €V and YorV 75, we have 
2.4.8 TO) = 5 te 2) Sn) 'g 2 (2n-N2P)4., (p) (x2) C4) 
= (any la’ (2-1) 2(p) Aes (p) Cy,-ayo"' 
+(4n)® a’? Qn) 2p) a, PY) a) om Naya 
-2§ (aea @w) '2(2N-/) z(p)o° (3-4)? 
c (4 5- a) f 2xop f= (20-1) 2 pean (a,-g)y F 
te | 3% (yap) -24& (n4,-h . 
~ Cen)"2 Cr-D2pA' yng F 
sep Een) e(pya’ 8m qa-a)y f. 


Having calculated the above, we are now in a position to 
choose the non-negative numbers & and Us Since 2.37 
is valid in the domain of 2.7, it is clear that WY 2o 


may be fixed so that 


2h9 TF ) <0 (x,y)é {x.y| ?Po?° va < $55 
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and 





2.50 J(F)<o (xsy)e { Xs¥ (P90 E 8 ¥Q-¥ € Shh 
By inspecting 2.8 it is clear that Le o may be chosen so 


that 
ge) T(F)<K<o (x,y) €} x,y |(x,y) € Wie ey as fa 


In order to establish that the auxiliary function 
is non-negative in V and is strictly positive in every 
closed set of V, we need only consider equation 2.39. 

Subject to the foregoing, it is clear that the 
theorem is established; however, as a final remark, we 


note that for certain differential operators we may choose 
2.52 GF =F (x,y)€ {x,y |(xy) € Vo<y-y ¢ 55} 


since T(F) satisfies 2.9 and 2.50. When this condition 
is applicable, one would of course choose © = 0 in 2.39. 
The conclusions of Theorem 2.2 shall be employed 
in the next chapter to obtain pointwise bounds for the 
solution of non-linear normally parabolic problems. [In 
the first section of Chapter IV the function % = Z(x,y) 


is not assumed to be strictly positive and 2.l\6 is replaced 


by 





el 


2.53 mane s 2 O 
V 
and 
2.5) Zi i ae (x,y) € D(y,) 


Under such an assumption it is obvious that Theorem 2.2 
is not valid since we would in general be unable to 
choose & so that 2.18 defines a strictly negative number. 
Because of this reason we use the fact that aij isa 
positive definite matrix and introduce a slight different 
auxiliary function. This new function is much more 
difficult to construct since, besides having a more com- 
plicated form, it is necessary to calculate a; (X09) 
for o<y<yp- Thus, even though the revised auxiliary 
function could be used in the linear normally parabolic 
case, the function as given by 2.39 would be easier to 
construct since the inversion of the matrix atJ is only 
necessary at the point péD(y,). 

When conditions 2.53 and 2.5) replace 2.16, we 


define A as follows: 


2.55 P= (y ae “t sup {-eNn-1) (p)5 em(ys-h) 





ee 
for (x,y)€ V where 
2.56 /p = as (Koay) (x =x5) (xIox5) 


For our present definition of F we notice that 2.35 is 


valid with 

2.57 p replacing p 

and 

2.58 Qing (x >y) replacing a. (p), 


In the revised form of 2.35 it is also necessary to add 


the following term on the right side of the equation 


ic 


2259 si %2(p) (2n-1) (8N (yg) 


“PY {- (2-1) zp) 5 (ancy-g) F | 


By exactly the same arguments as previously used, we deduce 
that there exist positive numbers 26 9, and Por such 
that 





é3 


2.60 J(#) <o (x y)€{x,y[5<2, 1,0 = Vo-vs 2851 
and 
rey He) <6 (x,yl€{ x. | PS 251208 Yor¥ § 254. 





Since 4(x,y) is strictly positive on the surface 
D(y,) it follows from the, conditions of the first section 
of this chapter that there exists a positive number 5 


such that 


oO N= =) 
2.62 53 gr) ( Dayrde- GM & <c 


We set 


B23 Ss = min | So; ce 


and notice that 


cl 


2.5), (L)<o (x, y)E (x,y | (x, EV, osyo-ys 28 § 


and 





2h 





2.65 I(L)<o (x, yEfxsy (x,y)EV, 0<€& y -ye28f 
where 

_ _ Mz! ‘ 
2.66 Leo (Ame) * anf {-Cow-) Bp) A; ANG a) f 


Since 2.62 is true, it is evident that 2.6 and 2.65 are 
valid for all 4? 0. By comparing Ce and X , we observe 


that we may chooseX =“ 4, such that 


207 T+ <2 (x5 ye{xsy \(x,yleV, 0£¥9-¥625f 
and 
2.68 THR) <o (xs WE fxsy |(x.9) €Vr0<es yy-¥s28,} 


Next we introduce the function @ given by 


fu 


2.69 Gr 


fy 


LS) - yp ay'h: 


"x40 L-(8N§,)'B (2w-1) 2(p) A$ 


4 


tae +z 
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where for the present os is an arbitrary non-negative 


number. To fix this arbitrary number we compute 
2.70 1) = (9'-Gy apt 
=) ; 
t (4 $)) 4 (AM-1) 2(p) a Ap G5 4-4) ses 


— (YN $0) 3 (en-i) =(4) a4q, Gy , 4) 
T (4M “ae a ; ait ( 2Nv-i)° 2(p) “Q ys 5 4) (> ) a x, 4) 


(nyt) 53 8 (aus yie2(en-) a(p) a8) Gens) 


<7 


+35 § a }~ (8v8,) 6 (en-) 2(p) AF 
Cad) = pl (BN 8.) 2 (2y-) Ap) A” $ 
+T(F+a). 


For $ % a $02 3. and (x,y) €V, the last term on the 
right of the above equation is a strictly negative 


number. Next we consider the following limit 
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2.71 Lim,  T(#F) 


Oo 


OS Asay S, 


i at _ (4a § ; 
‘1-(45,) 4 (2v-1) ZCp) +t 35 J 


(gg) 2 + Lm I(E4R), 


Po 


0% 4 £4,-§, 


By the remark following 2.70 and due to the form of 2.71, 
ded that th t iti 
we deduce tha ere exist positive numbers P o2 and eS 


such that for 


2.72 P =f, 


we have 


es I(F ) <o (x,yletxsy | A {9p 08¥S Fo 8 of. 


We are now in a position to define the desired 
auxiliary function. As was mentioned above, this function 
will be used when bounds for the solution of a degenerate 


problem are sought. The function is explicitly given by 





ey 


MY 


+a (4) € V OX ~4 SS, 







aa Por , OmeH> 8, 


+ 


ae = US) (gray SL gd) u-dep a BY: 


+ (4ns,d) 2n-N2(pe Ar +1 2 
“bp f- (ans) B, (20-1) 2¢p) 2° } 


4+GF tol fe =e . OS YS 3, 


On the surface D(yo)-{p} the function 4 ana LS 
derivatives are defined as the limits of these functions 
as (x,y) approaches (x,y,) with (x,y) EV ana XX. We 
now show that with a properly chosen the function ond 


satisfies the conclusions of the following theorem. 


Theorem 2.3. The function@ , defined by ies 
is twice piecewise continuously differentiable in x and 
piecewise continuously differentiable in y for (x,y)EV 


except at p = (X52Vo) and is such that 
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2.75 J(F) fo (x,y)EV 
Eno (A )<o (x, y)ELx,y| (x,y) € V,o<é¢ Yo-y$ 
and 
2.77 axi-$t = M 
vax {-¥] 


where M is a finite number. The coefficients of the 
differential operator J satisfy the requirements of 
section 2 and in addition the condition that the partial 
derivative of QD ng %o?¥) with respect to y is piecewise 


continuous on the line (x,,y)€V. The function % is such 


that 
eee min £ 7 o 
V 
and 
279 ZB(x,y)>0 (x,y) ED(y). 


Proof: The continuity and differentiability of & are 
obvious for (x, ye? x,y |(x,y) €V, o<yynys 3 of - For the 


lower part of V we calculate 





a 1F 3 25 -G>yS- 


Phy 
Oa OS a7 &o 


yp - (8/S,) (2-1) 2(p) yo t +e +o 


him \F} 


b Por 
P Po. 
O4 4 § Hom S, 
and 
- &, $= LC ay 
P=Pr 
06 $ Fo-b, 


UMS, 8 (2y-1)) Z(p) Qi» Cx 4) (x x) 


sept (enay'2 an-)eg ait +F%, 


Our task is next to select the number a 30 that 7 
satisfies 2.75 and 2.76. From the remarks associated 
with equations 2.55 to 2.73 it is clear that these in- 


equalities are satisfied on the set 
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2.82 LX 0h 





(x AE ve O< aot 8. ¢ U ix, a| on 5 oky edn Soh 
For fp >Po2 and o N57 o5 we compute 
2.83 T@)=1 2 Cy -afy'§ 1 Gav $,)' (2n-1) +16, 31- at, 692) ; 
ee 6.) 6 (5 Ae gerd een (E 
4 oN cs 3% S -*”) lay 34 zs. =x" ve 
¢ yee a4 2 
—(d -2) aoa 4,4) (* =) JO Oy) (x (4) oD } 
Lp \- Cens,y! a (2v-) 2p) a? } + {- (4-4) & t 
+g=' : (4-8) | 3 t4- (4NS, d) @n- “)E(py3 os 
+ (4ns,d) (2N-1) 2(p) (3 Ar a 3 


eyp {-(aus,y'B, (2m-D2part+ TE +a). 


In order to determine an appropriate choice of the 


constant li we note that 


ij ‘yo? S258) Sea er He 
2.8h a Asn (Xo9¥) (x x5) 85; (x, »V) (x8-x )7a ofl a> 
for all (x,y) belonging to the domain of definition of 
2.833; hence it is evident that a may be chosen so that 


2.83 defines a strictly negative number. 


With d so chosen it is obvious that the finite 





a 


number M given by 


2.85 = max{ -F } 


exists since Yr is bounded everywhere on V except at the 


point p and at this point one has 


2.86 ., qe 2 oO. 
i ae 


pev 


In future investigations bearing on the subject 


matter of this thesis, an attempt will be made to replace 


2.79 by 
27 Zoe es 


Next we shall establish the existence of an 
integral which will often appear in the bounds which are 


developed in the following chapters. 


Theorem 2.l). The function 4 is such that the 


integral 
2.88 zt Srey [Te] dv 
V 


exists. 
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Proof: As in the other theorems of this section, the proof 
is completed by discussing two cases. When 2.)\6 is satisfied, 
it follows from 2.35 that there exist positive numbers 35 
and (4 which are such that 
23 ie 1 
ae) eS 7) 2 (x ty).- 2 (0) (v9) 4(32N) (2N-1) 2(0)°? |. 
-1 
vexp{ F 8u(y,-y)] “ [(2n-2)2() 2°] } 


for ye ag and o< SaaS S54 éo: With ya replaced by 
Se it is evident that a similar result follows for Th in 


case 2.79 is applicable. Hence we have 
N+1 
a, —s— -1 2 
2.90 -d(%)? 5(7,-7) B(p)exp{ - lenty-¥) | [(2n-1)8(p) / } 


D < cy - ¢ ; 
for mele < (PRP and o te y< 6, 3. Next we consider 


= SET ®Y [Tp)]*dv 
= Stay | (3%) dv a Siyere)] [re,]* dv 


2.92 aes { (x59) | Psp £Py 2 gay < it. 


The first integral on the right of 2.91 clearly exists since 
its integrand is bounded over the domain of integration; 


also, since 


2.93 (8N) “th one2) <5 ; 





SD 


the second integral on the right of the aforesaid equation 


is seen to be finite by application of 2.90 and Theorem 2.1. 


3. Generalized Form of an Identity of Hormander 

The identity under consideration was derived by 
Hormander [16] for the purpose of obtaining a uniqueness 
theorem and estimates for normally hyperbolic partial 
differential equations. Essentially the same identity was 
employed by Payne and Weinberger [21 | to obtain solution 
bounds for second order elliptic equations. In what follows 
& slightly generalized form of the same identity is derived. 

Let ft(x,y), o cee: tN (x,y) be any set of functions 
which are piecewise continuously differentiable in x for 
(x,y)€ V. By direct differentiation and application of the 
divergence theorem, we have for any function ¥ which is 
continuously differentiable in x for (x,y) €V and whose 
eecona derivative in x is continuous in the interior of a 
finite number of subregions, the sum of which is D(y), where 
y is such that o<y<y, (in the above we suppose that D(y) 


is a domain of N dimensional space). 
ene uk nk L 
2-9 If lfoa *- fa, t- fa ‘ 4, Y5 Mp AS 
=-2 ff ery, ( (aity - AV 


a 


| 3 Ly n Ry 4 Ub 
+ SESE -f,4 fa +f “ag lK, % Y AV. 





The functions aij = alJd(x,y), igi = 1, 2... He Whateh 
appear in the last equation, are assumed to be symmetric, 
piecewise continuously differentiable for (x,y)€V and such 
that 2.7 is satisfied. (For the results of this section the 
piecewise continuity of the derivative of aid with respect 
to y is not needed; but, since this condition shall later be 
necessary, it is included here.) 


Let 


a 


WY _ ky 
2.95 cay LY 


Ps (5%) ES 


and then note that the vector 
7 ei y -| 
2.96 (T* 0) = (a Ty - Syl ne mg} I, 0) 


defined almost everywhere on $§ is orthogonal to (n, » ny) 
where in 2.96 we have made use of 2.5 and 2.7. From the 
above it is clear that (7J,o) is a tangent vector to the 
surface §, is perpendicular to the vector (0,1) and is 
defined almost everywhere on §. Following closely the 


procedure of section 2 of [21] we renormalize by setting 


i - wt 
2.97 (t,o) = (q* [sPIn in ]2 [a 7779 2, 0) 


where sy is the inverse matrix of aid, For (x,Y)E 0 


we have 


ied tee = 
2.98 as jt ty =a nny n 


where n is defined by this equation. Due to the above, the 


directional derivative on §S, 





35 
an zd i 
2-99 ye ~ Yt + Ty'0 = py ¢ 


J“ oe > 


Qy 


is in a tangential direction which is perpendicular to the 


Ny" direction. For (x,y)€ 5 we notice that 


ef hi “a Ly Sal] _y me 
2.100 a, 7 = A,o a " yv™ Ms Jum 
Mom st, 7 
=r Ly, av al 
wk 
“Ty 
R 
Also 
qigt_ tt 2 (Mee + hye) (aye 4% 
2-101 A 4! ; =a 4, %,— at aw )+ (38) ww 
_ kd 1 (24 
=Q Fe ae , 
Hence almost everywhere for (x,y)E§ 
by -72¥\* wWeyertre Fl -lipwy A 
Eea0c a Pag = M (35) + (47 Vee =m (3) +( | : 
rem \Vy My) -1l72 4% ( 24V44 4p | -i{pav\®_ (ay 4 
0,.0°7 Pode (SEY BIT Vo (38) 3) 


Moreover almost everywhere on 5 we have 
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aa a, f' tht =a wry sy Tr 


44 an 


ok wae 
=a, iT % ly av es 


k YY 
a ee + vw. Spy 


and thus 


sai oY : 
2.10) fim, DV ta, FORE =m toy. 


- 


Using 2.102 and 2.10) we rewrite 2.9 as 


| vy2 >Y 
2.105 JS ein, CSE) ~ (33 sy he 2a,ct't 3 33 341 ds 
fk yy 
~ . eSSS f Hp (4 Wald 


+ Sle, agi pha +P aly 4 dv, 


“R 


¥ 14 


We now choose c*, k =i, 2... WN; so that a bound for the 
square of the integral of the conormal type derivative 2.95 
may be obtained from 2.105. Previously we required that 

these auxiliary functions be piecewise continuously differ- 


entiable in x3 we now further assume that fin is bounded 


$9 
and has a positive minimum on 9°. Such a set of functions 
was considered by Payne and Weinberger in section 3 of [21]. 
From this reference we note that, if 5 is starshaped in x 
with respect to some line perpendicular to the x plane, then 


we can choose this line as the x origin and take pk = xX, 





Bi 


Also if V is a domain between two boundaries each of which 

is starshaped in x with respect to the above perpendicular 
line and if V contains a cylindrical shell of radius re 

from this line then we may take fk = xk (rer)- For further 
information bearing on this choice of the auxiliary functions, 
the reader is referred to Appendix A of [17 ].* With gt as 


defined above, set 


_ eo 1 
2.106 Mm = saa n £m, 
Since aij is positive definite, there exists a constant 


Cy which satisfies 


c ie. ee hale ee 4 
2.107 ~{f,ya ~f,, a “f,78 + f a, en by ts 4 
fer well (xyy) ¢V- (C, igs an upper bound for the largest 
eigenvalue 6, of the coefficient matrix on the left with 
respect to aij.) For a crude upper bound for eS we define 


2.108 gid = -|t,ka*J-2, asker, Jat : ta, | 


then 





666 also, Bramble, J. and Hubbard, B., Some Higher 
Order Integral Identities with Application to Bounding 
Techniques, in print. 
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2.109 ex sma Bop - z _&, 
Dove TES FOES 
ihevyeepys 
a,>_ §. 
7 A 
B 
<a = 6 
where 
5 =n a 
Zed. 0 B max [> ) et 


Substitution of 2.106 and 2.107 into 2.105 yields after 
applying the indicated inequality to the resulting expression 


2.111 (i-¢) m, S(3 Sy) ds & aS (mF rng) (ty fe at) Ls 
Sof in( BEyids +2 SSS RY (alty, 
¢ A y A a?) 
+ C SS ay, 4 dv 
where & is an arbitrary positive number. By setting 


Peeinie c= 


we may write 
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ams S(SbYds «H+ h, SSS He, (aMy,,),, dv 
on JNJ a4 ey ¥, AV 
where 
2.11) i =a. S Ont ing) (Qe, fee SS) ds 
\A 


in, Sfmt! fom, ( BE) as 


and clearly 
2.115 AL = QO 


pag ty Y ="ortor Gwwy)eS . In the following chapters great use 


shall be made of 2.113 when the boundary value problem under 
consideration is of the Dirichlet type. 





CHAPTER III 
NON-LINEAR NORMALLY PARABOLIC PROBLEMS 


1. The Non-Linear Mixed Problem 

In this section, bounds for the solution function of 
a‘boundary value problem, which is normally parabolic, will 
be derived where the solution function is initially given 
and thereafter satisfies a non-linear Neumann condition on 
the boundary and a given differential equation in the 
interior. 

We assume that V is an N + 1 dimensional domain with 
boundary D(y,)U SUD(o) which is such that the conditions of 
the first section of Chapter II are satisfied. Our goal is 


to determine pointwise bounds for the solution of 


J (w) = (a*w,,),.- z [ = F (x,m45 Ww, vw) (x) EV 

— CO W VW CoO 

ee W(x,0) = 49 («4 ED(2) 
W_ 4 

85 = a "wom, = L (4,7) (@&A)ES 


—O < W < 0 


where W is continuous in VU D(y,), continuously differentiable 
in x for (x,y) €V and is such that its second derivative in 
x is continuous in the interior of a finite number of 


subregions, the sum of which is Diy); where y satisfies 


ho 





Wi 


Ey) SIR : Also we suppose that the first derivative of w in 
y is continuous in the interior of a finite number of 
subintervals, the sum of which is the line (x,y) for every 
fixed x where (x,y)€V. In addition it is necessary that the 
solution of 3.1 satisfy 

lim w(x,y) = g(x) 


(x,y) (x50) 
(oy) EV 


i ~ = 
lim aid (Xsy)ws, (Xsy) = aid (Xsy)Ws , (Xsy) 

(x,y) —> (apy ) 

C35) SV 4=/,R,....N, 

(xsy) ES 

The boundary value problem is assumed to be such that 

the components of the symmetric matrix aij = alj (x,y) are 
piecewise continuously differentiable in x and are such , 
that 2.7 is satisfied in V. In the closed set V the function 


Z= (x,y) is piecewise continuously differentiable in y and 


satisfies 


> Oo 


3.3 min Z(x,y) =n, 
V 


The boundary data f 9 L and 4 are assumed to be integrable 
and square integrable over their domains of definition for 
any bounded continuous function w which is continuously 
aifferentiable in x for (x,y) €V. In addition f and 2 
are required to satisfy a Lipschitz condition in all but 


their first two arguments. Hence there exists positive 


~~ 


numbers M,, M,, Me, i = 1, 2... N such that 





2 


~~ ~~ 
MAR ty w, ,ww2)1§ M, lw,-w,l+ Jw eM | 
| Gm)eV 5 co CW, VW < CO 


3.44 


Reg wi) Lea Ws) S M, |w,- wa | 


(Kj WES, —o<wK< oo, 
= 
Except for the fact that L = 2 (x,y, w) is required 
to satisfy a Lipschitz dondition in w, the problem of this 
section is a generalization of a problem whicn Friedman 
considered in [11]. In that paper the elliptic part of the 
differential equation was the usual Laplace operator, f 
was identically zero and Q was a monotone non-increasing 
function of w. Under these restrictions the solution of 
the boundary value problem was shown to exist by means of 
the Schauder Fixed Point Theorem [26]. In a concluding 
remark the author indicated that in a later paper he would 
attempt to extend the existence proof to general second 
order parabolic equations. Whether or not such an extension 
has been carried out is not known. The existence of a 
solution for a similar problem was however established by 
Piskorek [22 | and also by Pogorzelski [23]. 
Let p Ebi.) be the point at which a bound for the 
solution is desired. We assume that the functions a‘J ana 


% satisfy a Lipschitz condition at p éD(y,) for (x,y)€V; 





that is, we assume that 2.8 and 2.9 are satisfied. In 
addition it is necessary to suppose that f(x,y,w, Vw) is 
bounded in some neighborhood of the point p which is con- 
tained in V for any continuous function W which is continu- 
ously differentiable in x for (x,y)€V. Hence for any such 


function w there exists numbers § 


>? [53 and ¥ (w) ali 


greater than zero such that 
3.5 f(x,y,w, vw) < /6(w) 


for f? <P> and a 6559 <7: 

The desired bound for the solution is obtained by 
choosing an arbitrary function Y (x,y) which is such that 
Gi) eC and YW (x,o) approximate f(x,y, ¥,7 ¥), Lizsys ¥) 
and 4 (x) respectively where we further require that D be 
twice piecewise continuously differentiable in x and 
piecewise continuously differentiable in y for (x,y) €V. 


Let 
346 Tey) =w(x,y) - (x;y) 


and compute 


1G. Fx 4, w, Vw) - Pes 7) + Fee) 


(4) EV OX WOW, VP <00 
3.7 VY (0) = GY 65%) € D(o) 
35 = X Oy) L (> (x40 u) +L& 4) 


where (Ky) ES 5 ~o Cw, poo 





3.8 Pix,y) = F(s;y, Pa vo) — 3 7) 


and 


3 
329 L(x,y) =X(x,y, f) - Ae . 


Due to 3.3 and 3.11 we may write 


3.10 CY) <i Yl eATY | + FG) 
and 

y 
eel u< M. IY] + Lb Om). 


With T(Y) ana Yass defined by 2.6 and 2.10 respectively, 


from the divergence theorem it follows that 


3.12 i LY, T(Y)- ¥T(8.) flv 


= ff {21 35 Sl y BP} ds = on Sp Yels 


S(*4) 


~ SS ang zvds + Ll Sp BY ats 


Sis) 4 
where 
3.13 Mie = { xs) | (sm) EDC Ds o<m < yf 
and 


3.1 Sy) = { (| (99 VE Bly V-D(4YIs 0 CT}. 
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In 3.1) we let y approach Y. and use Theorem 2.1 to obtain 


° 5 xK = , i _ 
3.1 YK, Mo) Bors Lys {YT (p)-% T(Y) } dv 
+ StH, BY BF bs - Sn a2 vd Sf 2 Vas 
which becomes after substitution of equations 3.7 to 3.11 
3.16 Iwi.) - Xp) 
7 WN oF (ey )aV — a WN. f(s, 4,4, ow) aa 4 wr otdV 
t iS ue LL (x, yy w)- RO, es dS 
- SF 54S ~ SSm yp 24 ds | 
<[dlny SLY F Cap) dy [+ ae Nyt foramen) toyed dv| 


+4 SY ‘dsy* KR ser “ zy’ ds +M MS a2d sy) 


where 
3.17 X(p) = ~O%,4,) + JS LOx,) Up ds 
+ LU Mpecds - for SSS Xo F ay. 


V(4) 
The terms on the right of 3.17 are all known; 


however, it is necessary to prove that the last integral 


converges. For this purpose we write 
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O°Go VY) 


=i 1 Sf YpFdv + S&F dv) 
G7'te | W§)-V(fas8,) V4) AV(A,8,) 


== Sf Spt (Ku, TOV = ina %pI(e) dv 
v-¥(2.,5,) v-ViA,§,) 


t Lim SS %p{Fe 4,74) -I(P}dv 
$770 VinOV(A, S.) ‘ (Os 


where 


3-19 V(A 5 51) = iesd| es ? osy.-4 <5, f 


and on the surface D(y,)-{p} it is clear that ¥ » is defined 
as indicated in Section 2 of Chapter II. We now consider 
the three integrals on the extreme right of 3.18. The 
first exists since J D is bounded in the domain of inte- 
gration and f is integrable over V. The second integral 
exists since its integrand is bounded over the domain of 
integration. For the last term we use the fact that 
[£(xsys Y,VY)-JI(¥# ) | is bounded over V( Pos 35) due to 
3.5 and the definition of J; thus the integral exists since 
Y, is easily seen to be positive and integrable over 
V( Pos B5)- 

From the above remarks it is clear that all the terms 


on the left of 3.16 except w(x,,y,) are known and hence to 
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obtain a bound for the solution it is only necessary to 
estimate the expressions on the right. To obtain this 


result, we introduce the function Y defined by 
3.20 Y = ¥ exp [-K(y,-y)] 


and calculate 


3.21 I(¥) = 3( ¥ exp [-K(y,-y)]) = 5(%) exp [-K(y,-y) 
where 
Bae2 3(V) = 3 %)-k2? 


Let ri, i1=-1, 2. ..N;, be the set of auxiliary functions 
which was introduced in Section of Chapter II, that is, 
the ri are piecewise continuously differentiable in x for 
(x,y)€ V and are such that 
Baas o<m, = min ari. 

g 


With the ft as thus defined we may write 
suk | aR a 
3.2h YJ ¥ds < omy Sg f aAabes 
5 
~ & mw { n~k& iC 
a) SSS # yt dv + om SSS Y fis dv 


«Mi, ae 
: omy 4 Salty 





A 


D a na) ~R 
vi dvt Com, + ae TY av 
where % is any positive number and 


N 
3.25 M, = max | > ete! } 





18 


Es at 
3.26 M5 = _ {2} | 


In arriving at 3.2lh we have made use of 


ne | 

3.27 (fy ys Bt, . 
which is easily established by expanding 
3.28 > (#'y at 


pa 


From the divergence theorem it is clear that 


3.29 Sif ¥ T(®dy = “SSS 4, a "G.dv + SSG ay 


“Sif red “dv “af “zal s refer weal 
D(4 9 


re OH) dv in 2 Os. 


Using 3-/ amd 3.20 to write 


3.30 F(P)< FH + ANE + Fg) soup |K(yo)] 
and 

yy q N “ = 
3.31 ee el + LO 4) ope [Keg] , 


we have by means of the above 





Lg 








3.32 SS 9 F € SSSIFH T@ dv 
cH, SSO y + et SiSatty,, 0 dv +e, LP 
0 Vv 


+ SSS Pdy 4 Be, SUJ F* anf [Rk (yy-4)] dV 


and 
as, a 
e SP Sy dsl ASP ds + BSS P'ds 


+ a, SS LE amp fancy-g)) ds 


% 


where Los o., and a are arbitrary positive numbers and 


Ww 


NV 
3.34 M, = > crty?. 


Substitution of 3.32 and 3.33 into 3.29 yields after the 


necessary transposition 
Ma ‘bo * mye 
3.35 Ch Ra) SfJa'? &, tv + hong My =F Fe - PMY? ‘AV 
(M+ s+ zm) SS? ‘ds +4 zie 2 up any, As 


+e, SIS Feup Rig g)|dv+ xq SComlengealds- aS Ved 
where (4 


3.36 m, = max 1 34 


and 
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3.37 Me = = {-n a}. 


Except for the first and last integrals, the right side of 
3.35 is completely known; the last term may, however, be 
neglected since 4 satisfies 3.3. In addition we note that 

the left side of 3.35 has the same form as the right side of 
3.23; and, after a close scrutinization of these inequalities, 
it is evident that they yield a bound for the integral of the 
square of p over S provided the arbitrary positive constants 
Cs and K are properly chosen. In a particular problem the 
best choice of these Sepiipery constants would depend on the 
functions F, G, and L and the M's appearing above. For 
definiteness we assume MF © and make the following selections. 


(If M, = 0, obvious modifications are made to obtain the same 


B 
result). Let 


3.38 of, = i, 

3.39 a = 1 

3.0 = gt Me + 3+ BM) 
Bend a, =1 

and 

3.2 K ra : + My + + . 


i. 2(2Mo+14M-)M, 
+(2M,+14M,) a + 2 ; = K 


2 oe, mh 





Sl 
where we have assumed that 
= 1 1 


If this is not the case, then one need only choose a SO 


that 
6 
i 
3h M., +t +d M,> 0 


and the same result is now obtained by appropriately defining 
K and A,. Substitution of the above into 3.35 and 3.2h 


gives 


5 b+ds < (2%,+4 Fi ‘ 
3-4 SSPds < QM +14!) Lf ¢ 2 Lyp (2K 4.) abs 
+ SSS FS sup [ak G.-4) Jel v SS Loe fan, aldst 


where K, is defined by 3.42 (assuming 3.43 is valid). We 
use the last inequality and 3.20 to obtain if K,<0 


3.46 SPds< GAH AY Lifer ols 
5 (o 
+ SJ Fae 2-4) dy -} SS Pegg -2k 4) dst 


and if K,;20 then 


3.47 SS Fads (241+) NG (axr4,)als 


- SSY F exp 2k, Yor d)\dv + {Ceyphen,-a)Hs §. 
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The right side of the applicable inequality is known since 
F, L and G are square intesrable over their domains of 


definition; hence we may write 
3.1.8 J pls < 


where $b is a known number depending on the function Y. 
After having calculated 3.1\8, we ncte that a bound has been 
obtained for the last term on the right of 3.16. 

If a function I|'p can bs determined which is such 
that the conclusions of Fheorem 2.1 are valid and in addition 


the condition that 
3.h9 Fp) = 0, 


then it is clear that |p is a fundamental solution of J and 
the first term on the right of 3.16 is identically zero. 
Dressel [3] and [| established the existence of such a 
solution when the coefficients of the differential operator 
are sufficiently smooth. From these references it is clear 
that even when its existence is guaranteed such a function 
may be quite difficult to construct; hence we shall assume 
that |'p is not known and that the first term on the right 
of 3.16 must be bounded. To obtain the desired bound we 
introduce the auxiliary function Tr aerined by 2.39. Since 
the differential operator of this section is identical to 
that of Theorem 2.2, we conclude that oO and & may be 


fixed so that 





Be 


3.50 Ge 2 o (x,y)EV 
3.51 J (%)<0 (ayy) EV 
and 


(F <0 GF >o 


) 


C4 


Babe 


for every €70 and y satisfying O<€< Toy Soi oe Next we 


define the adjoint F of J by writing 
3653 J(%) =T(¥%) = xe 


Because the singularity of % atop is of a slightly 
lower order than the singularity of the parametrix Son: it 


follows by the proof of Theorem 2.1 that 


3.5) a a CL SZ ids 
47> Dg) 
4 <4o 


= 10) 


since 4 is strictly positive and & is a non-negative number. 


By means of Green's Theorem we write 


3.55 Dre SISA F(P)- VF (F} dy 


na 
: Say pd - Singha Ps 





5k 


where we have used 3.5l and the fact that 


3.56 


for pé Disa) : 


Lom SS Je 3S - ot 2, 
424, SU) 
4 <i, 

= SSi% a0 pe = tds 


If we assume that K is non-negative, then 


expansion of 3.55 yields 


3.57 


Doe iSSie F(t) - VI) 
a 


+244 vat yi, +h her pe} dv 
Sf DF)" eng [2k (y,-a)]|ds + Zi SS Pads 
taM, SH bd: - Se oF ds = SSnG 2 tds 
+f b G2 mp Qk.) ds 

¢ SSCFLY snphan gays + MTS 2np (2ky,) 


— Zep (aky,) dS = y 


where in the right inequality we have for definiteness set 


o = 1 and 
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as ~oe 
3.58 Me = = 10, 1+2M ot ~ am aia gh. 
By means of 3.7 and 3.21 we have 


3.59 Lhe TCP) Av = RANG P Fon [KC4,-4)] dv 

eI OFE {fey,.w, 9) ~F64,0, 70} axplitlana)] AV, 

4 
Also 
3.60 2 SSf HIF ae dow, 7) 
v(4) , 
-f& 44 v €) | 2p | C4o- )] dv 
S25 SUF Md 2A FEE dv 
V(-4) v4) 


(aR, +4) FEN + Be (Cale e G. dy 
V(4) a = 


where M., is given by 3.3 and a is an arbitrary positive 


number. [Let 


M 
2 Mu Boe: 
3 61 5 2a, 


then by choosing K so that 


be 


a t. BD 


O 


1 ne 
3.62 ea a | 2ii, + 


MN 


and using 3.60, it is evident that 3.57 may be rewritten as 


3.63 Dim f Ste Pinply-a-PIasdv ¢y 


37do Vig 
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With “f defined on D(y,) - {p} as is indicated in Chapter II 
¢: a 

and with the constants O 5 and ls chosen so that ia fs and 

5 a8 &>> we have from 2.35 for appropriately defined constants 


Css i = iO’. ° @ « L3> 


at Lom SL REV F aypo [ig (m,-g)|- 9° F 
‘oh {iy i PYF cup (g-d|-P TDS av 


= ie §fs 12h |= 
47> V(A,2,)- Vey (x, 4,0, 7) 


- T(v)| 2xp[ k, G4,-41 ~ ps T(F)} Ry 
~ do +) 1 cia a a ( o” Te Co ia ie c 
17% V(P,,8,.)-V(4) ‘4 4) in SG 4) 2 


Eis (Ay- 4) =| dy40 (- ea, M (.-74)) | len-) zp] Y} dv 


In arriving at the last result we have integrated by parts 
first with respect to 7% and then with respect to y (see the 
proof of part a of Theorem 2.1). From the above it follows 
that 


3.65 [S12 % Fawp [yd ]-0° F (Fd v < \, 


or 
From Schwarz's inequality and the fact that 7 satisfies 3.51 


it follows that 
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3.66 - LSE FF ou {K, (45-4) A V| 
> {SS OIF) dv}* 


LSE Iq Rte fee, c4-0)] dv}? 


The factors on the right of 3.66 each exist. The first 
exists since Y is continuous on V and is hence bounded on V 
while T(F#) is easily seen to be integrable over V. For the 
second factor we note that the singularity of 4 : at pe€ D(y,) 
is less than that of J( ¥p) 2 and the result may thus be 
obtained directly from Theorem 2.1. Due to the above 


3.61 {{lJ-3 08) Pov ES | SES inp fing ol s* 
{= SIS FTG) dy 4 ELIF P Fang I, yo--4) Jal 
SE IF" FF) oop fey ,-d] dv 3 


and then, after substitution of 3.65 and 3.20, we write 


3.68 {-S{6 TF) dv} 
LILO] @ F)* anole K, (4-8) dv} 


19, + SSEIGY ry ope y,-a)dv¥t 
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Because of the form of oa and Y ps we use Schwarz's 


inequality to obtain 
3.69 Dim SSS ¥ F(¥p) olv | 
a 00 ¥(4) 


2 


< {SSL I Toy 
SSSI T A cap Beaty a) v)* 
19+ SSE) OP Yioep [tm c4- gard 


where the first factor on the left exists due to Theorem 2.) 
and the second factor by the remark following inequality 
3.66. It is clear that 3.69 is the desired bound for the 


first term on the right of 3.16. 
To obtain the desired bound it is now only necessary 


to bound the second term on the right of 3.16. We write 


3470 Bon i Sp {fw ew)-F (x4 Y VYYf lv | 
< Lun (Sf &, } M lvl + ily, | aly 
$°%G> Vg) 
<iNy Xp F open ged Avt Uae SSH Feat akan Wilavt™ 


LSP owpakengea yd ‘Le Spa anys 7 <'e\ag 
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In the last inequality the first factor on the right exists 


since for ah appropriately defined 


ar SIS F axgo Ak(y-a)jov < Cy + 


"6s ‘ if (4, -4) ® “7 2p [BNA (4 Dl | (ew- newalta daduody 


for & sufficiently small. The integral in the above 
expression is seen to exist by means of integration by 
parts as in the proof of part (b) of Theorem 2.1. The 
last factor on the right of 3.70 is now bounded by recon- 
sidering 3.57 and 3.59. From these inequalities we obtain 


3.72 Y 2 Ln ff ( LF PI(¥)-% IF) tPF Gary +2K2F ty 
“fo VC) = 
~<a, 
On Sf t- FF enp[e K¢,- 4)| 


470 «W(4) 
4 <Go 


re (Zk 2 =o, - 2H, ~0,)% F(z + a fray, 5 dy. 


In the last equation the choice of the ot. depend on the 
magnitude of Ma.” and the boundary data. For definiteness 


we set 
3.73 Oo =1 
and 
1 1“ ~1 
a7 (oem o, = 2M3 (a, -M4) 7 for 5> 3% 
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or 
a =] -l1. 1 

: b = = 
os 3 a, M2, for Me > 5- 
Next we choose 
Bef) We K = i (ascit a +5M, (a, =M 3) “teati,2 } = K 
or 

ages 2 
Bei D K = sa {142 0 te M, “it } K, 
whence either , 
3.75 a Bem SSS 2A FR + 2ma'Gaty & A 
“so V(-4) “ 


<), + SUF E weplanty-pldv = 5 


Ee ieeaioy ig be me 
3.75 bes Slam Pe OMe ba 9.8 \dv 


< (21,45) 19), SFE expe ks (go-—)|4 V8 = 


where is given by 3.57 with K = K,. 

Subject to the remark following 3.71 it is clear 
that equation 3.70 and the proper form of 3.75 provide a 
bound for the second term on the right of 3.16. 


By combining the results of this section, the solution 


of 3.1 is bounded at the point pED(y,) by 
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3.75.5 I worse) - XCp)| 
<{ S{sEr FY tev? 
= CTA) (FF) aep [2 ky (,-¢)} dv $* 


+19, t Sifts "G FF) up [ek, (4,4) jay 
LY F ep Lae lq-adav § “f ) , 


Bi QU Gz ng ds} +AU S Ads *) | 


In the last inequality the functions ¥ and “ are given by 
2.10 and 2.39, and the constants K, by 3.2, 30962, aad 3.7. 
The computable bounds oy Ne and 9 are given respectively 
by 3-48, 3.57, and 3.75 while the function X(p) is defined 
by 317. 

If problem 3.1 admits a solution w which is twice 
piecewise continuously differentiable in y for (x,y)€V, 
then theoretically the pointwise bound may be made 
arbitrarily small by improvement of the approximation 
function ). In this regard we note that the right side 
of 3.75.5 is zero if F =L =G=0. Furthermore, due to 
the form of the bound, it is clear that the result may be 
improved by means of the Rayleigh-Ritz technique provided 


the problem is linear. 
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By referring to the proofs of the appropriate theorems 
of Chapter II and the methods of this section, one observes 
that the boundary value problem given in 3.1 may be 
generalized by allowing the components aij = aiJj(x,y) to 
be non-symmetric. Such a problem was not considered since, 
even though the techniques are exactly the same, the 
individual inequality expressions are naturally much more 
involved. In the next section, where the matrix alj is 
required to be symmetric, a method is given for synmetrizing 
a general problem. f 

As a second generalization it 1s easy to see that the 
differentiability of the alJ, 1,j =1, 2... N may be 
reduced by merely requiring that ayy and a, 4 be continuous 
in the interior of a finite number of regions the sum of 
which is D(y) for o<y<y,- The atj are themselves con- 
tinuous in D(y) for o<y<y, and piecewise continuous in V. 


Under this reduced differentiability, one may easily show 


that the techniques of this section are valid. 


2. The Non-Linear Dirichlet Problem 

Our aim is to bound the solution of a normally 
parabolic problem which is such that the solution is 
initially given and thereafter satisfies a non-linear 
normally parabolic differential equation in the interior 
while assuming given values of the boundary. 

We assume that V is an N + 1 dimensional domain with 


boundary Dive y 6 0 Dioy which is such that 2.) and 2.5 are 


“¥ 
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satisfied. The boundary value problem may be stated 
explicitly as follows: 
s(w) = (ati w,,),, - eater (x,y,w, 0M) (x,y) EV 
| J J 00 <W, TW <9 
3.76 w(x,0) = g(x) (x,y) € D(o) 
h(x,y) (xsy)E 9 


If 


w(x,y) 


where w is continuous in VUD(y,), continuously differentiable 
in VUS and is such that its second derivative in x and 
mixed derivative in x and y are continuous in the interior 
of a finite number of subregions, the sum of which is D(y) 
for every y satisfying OS FSG: In addition it is necessary 
to assume that | 

lim w(x,y) = g(x) 

(x,y) (05x) 

(x,y) EV 


eo 
lim w,;(x,y) = Ss (x) {1 =a 2 . . eee: 


(x,y) (05x 

(x,y) EV 
The differentiability conditions which the solution of the 
present problem is to satisfy are obviously much stronger 
than those which were required for the solution of the 
boundary value problem considered in the first section of 
this chapter. Additional conditions must also be placed in 
the coefficients of the differential operator; in particular, 
the components of the symmetric matrix aij = atJj (x,y) are 
piecewise continuously differentiable in V and satisfy 2.7 
while the function % = 4(x,y) is piecewise continuously 
aifferentiable in y for (x,y) €V and satisfies 





ou, 
2070 min Z(x,y) = Ms >0- 
V 


The boundary data f, g, and Bos are assumed to be integrable 
and square integrable over their respective domains of 
definition for any bounded continuous function w which is 
continuously differentiable in x for (x,y)€V. Moreover, 
the function f satisfies a Lipschitz condition in all but 
its first two arguments. Hence there exist positive numbers 


M and Wi, i =1, 2... N, such that 


: Cissy wae 7 we ) fbx w4,7W <M. |lw.-w +i | w =-W 
B79 «= ayy, ys VW )-f{X,YsWos 2)|< | fl 2| | Nes fon 
for (x,y) €V,=0<w,7x << ©, 


We remark that the problem defined above and the additional 
problem which we shail obtain from it by introducing an 
arbitrary function are very similar to a problem for which 
Friedman, [10 ] and [12] » and Zeragiya, [27] and (28| 
established the existence of a solution. 

The desired bound for the solution functicn will be 
obtained at the point péD(y,)- As in the previous section 
we assume that % satisfies 2.9 at p (notice that 2.8 is 
satisfied due to the differentiability of the atJ) and that 
f is bounded in some neighborhood of p which is contained 
in V for any continuous function w which is continuously 
differentiable in x for (x,y)€V, that is, we suppose that 
3.5 is satisfied. To obtain the desired result we choose 
an arbitrary function WY = (x,y) which is such that J(/)) 


and W (x,y) approximate f(x,sy,Y,7Y), g(x) and h(x,y) 
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respectively. In order that the divergence theorem may be 
applied in what follows we further require that Y¥ be twice 
piecewise continuously differentiable in x and piecewise 
continuously differentiable in y for (x,y)€V. Also we 
require that the Vs, be piecewise continuously differentiable 
in y for (x,y) €éV. 

Following the techniques of the first section of this 


chapter we define 


3.80 Y (x,y) = w(x,y)- ¥(x,y) 


~ 
= 


and then compute 


J(Y )=f(x,y,w, gw)-f (x,y, , TH)+F (x,y) (x,y) EV 


3.81 Y (x,0)=g(x)- Y(x,0) = G(x) (x,y) € D(o) 
Y (x,y) =h(x,y)- Y(x,y) = H(x,y) (x,y)é 8 

where 

3.82 F (x,y) rs P(F.Gak ; VY)-T(Y). 


For later use we note that 
3.83 ocY) |< iM [4] sit 1 y. | + be(x,y)| - 


With J(Y) and Cup defined respectively by 2.6 and 2.10 
we note that 3.15 is valid for the present problem. Sub- 
stitution of 3.81, 3.82, and 3.83 into that equation then 


yields after transposition 
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3. 8h | WGw,4,) - Xp) 8 ) Av 


4 <4 


He sip Yp [£& 4,wew)- thy, ev} dv 
47qo 


‘t § fo 


HIS MPL Sag as} 


where 


3.85 K(p) = aes) JS H Sds be «NS Fy 
‘to Vly 
fc <4 


+ SS ¥, 2Gds - ae y¥ 2H ds. 
4a 


D(o) 


The volume integral on the right of 3.85 is proved to 
exist by the argument associated with equations 3.18 and 
3.19. Hence w(x 5sy,) is the only unknown term of the 
left side of 3.8) and our desired error bound is then 
clearly available as soon as the expressions on the right 
of 3.8) are bounded. Since the unknown volume integrals 
will be bounded in terms of known quantities and the 
integral of the square of the conormal derivative of ¥ 
over S, the last term on the right of 3.8) is considered 
first. For this purpose use is made of 3.81 to rewrite 


2coll 3eag 
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3.86 SS (33) ds 
a Sf an YAY 
+ in, Sipe f VY [2 34 a a FO 4, w, 7) fey 7) )+ Fos) Pel 


where MW ts given by 2.11. Since f satisfies 3.79 we may 


wite 

3.87 fis" ee + F(x a, WIW)-FKs 4, 78) + Eloy) al 
<2, Sf “(PMY av + Bl (si )dv +4 LI (Py palv 
+h LN Vay ba Ye ldy +h (Ae Yl 
+E SSS OY dv + ZS FF dv 
< (rae, + 3 + wa) Say, Ye dv + Sf v"dv 
+e SSS (Sgyav + + SS Fd 

Suse 


2 — 2 
3.88 M Bn \3 2_(f4) } 
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and a, is any positive number. In the last inequality the 
constants M; and M3; are given by 3.25 and 3.34 respectively. 


Hence 3.86 becomes 


3.89 (Ss ods cA (+b + oe we) Las, dv 


BEE SIS 9 ait j dv + & Lite lV. 


Due to the differentiability of the solution of 3.81, we 


may write by means of the divergence theorem 


3.90 {Sf . T(¥)dv <—SSf =" at t dv 
V V 
+ SS 3 at ds - on, JSS (35 )* dv. 


Application of 3.83 to the left side of the last inequality 


yields 


3.91 eS ridv ls SSS sgh fel A lyfe lelt dv 
BSS Ya +e SN vay + BIB Ydv 


V 


See 


+ Zea, Shay, b dv t Bi (Sa)dve 1 LIS e* ov 


where the ef are arbitrary positive numbers. Subsequently, 
when we use 3.91 to rewrite 3.90, the arbitrary constants 


will for definiteness be defined as follows: 
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m. 
a= — a= co 
ae92 a = - = t = - 


The first term on the right of 3.90 is bounded by noting that 


3.93 - I § yt alt y te dv=-k Om Sf ya'ty, ds 


424, D4) "4 


tz SJ G,, Q af ds + ae SST 4 1a by, dv 


D(o) 





meek 
x SSm a" yt obs 
where 
N 1 
is dats jagils)s3 
3. 9 ™ - nl 3a 


Using 2.102 the last term on the right of 3.93 may be 


rewritten and bounded by 


3.95 ~ Sn a'4, tds =-Z Soin (S5)ids-b Slog (38a 


cB SW ds - 4 eon, Eds 


where 


3.96 My = max {=n tng} 


To estimate the second term on the right of 3.91 we first 


write for an arbitrary positive number As 
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ale 
Ae 


son (f BMds BIS Bids ck, SRY A 


S 


and then, almost everywhere on §, define the vector 


3.98 ~ (pt, , . . tt, a) 


where T is tangent to the surface 8S, is orthogonal to the 


vector (t2,0), given by 2.97, and has unit Euclidean length, 


LoGe 

3.99 s (tty? 4 ov 2 

We denote a directional derivative in the T direction by 
3.100 st > y T+? ay Tt 


Since the space is Euclidean, 


3-101 nt = ns for every i 
whence 
N 
UY 2Y LE Dy 
= =~ t\R Se 
BpLoe ON, mal 2 Ym, Tea on s (- me) oT . 


Transposing and squaring 3,102 yields 





TL 


Aes (1-2 )*( a =n 


NV 
Sy 5 oY 
hy ay t Cm) oT 


< AM (> m2) aay He t ae (35) 
= hn (S mz) ( am) (3 ae t a (> *)( an) (32) 
+ a( ray (SEY 


Since 2.5 is satisfied on S, 
2 


and hence from 3.103 


3.108 SS ds <n, (SS IAs + fe Leml (Bids 


+ ise om [a,m(i-mz)] aes 


where 
106 1. = max{2n° [2 n(1-n°)| “t 
De 10 — g yloo y ‘ 


By applying the above results to 3.90 we have 
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3-107 SIS HY V cp JIS ¥%dv 
(di, tie) SIS a" w, dv 
s / 
+ ( (Zz? + z+ wa) SS ( of )'ds + im, Sf Fy 
is Sy G,, ate. ds ~ HSS tony (2 )* 4 
t x Mme ye m Cray *)) (32) ds 


Zz jf Cm) (3+ ds ~* Ose Saag ds. 
p ONG 
<p 


In 3.107, let 


which minimizes the coefficient of the integral of the 
square of the conormal derivative. For what follows we 
assume that (Mg +2VM, 5) is a strictly positive number; 
when this is not the case, obvious modifications are made 


to obtain the same result. We set 


a ee -1 
3.109 a, = Te (Mo + 2 VM, 4) 


and then rewrite 3.89 as 





13 
3-110 I (BiYads < 2H my (6, + Cam, (1+ 20M) 
rns at BM, eV) (dees, +28) } Slay 
ef ae + £2 (4tevm, tS Vd + om ie 
t2 (M+2VM,) tear mm, SiS dv + ELS, a ‘a ds 
4S min (Se)ds + Hm, Sn lam engl (QE) ds 


va, SS C-may (35 Jas Lam [ma Beds $, 
ee 


if the Last term on the right of 3.110 is neglected, then 
the only unknown terms on the right of this inequality are 
the second and third. To estimate these unknown terms we 
make the transformation of variables given by 3.20 and 
observe that 3.29 is valid for our present problem. In 
place of 3.35 we write 


{yO AW 
pane (1 BE) SSS ae a 


+ (km, - My > M a, 3) JJ dv 


; #8) 1 om Yop lingrdps-2 I Ve dS 





an 


Po G 2 ou40 (2ng,) ds t Ra, NY F oxy Bage-g) lV 


D(») 
me {Sm2 Ho pk (4,-4)]ds 


In the last inequality, let 


a 


© 
3.112 q. ~ i. 
3 
14 AL, =1 
35.183 3 
M | .~™m 4M 
3-114 K = m2 > tM i+ 2a, ee Y F |e, + dzjmom_(My + AVA) 


S m™- -| ONG 
‘4 +e mr" (4, th Vo) rae + Fe )) Es + 
3M, a 
* Zang(My +2 VA) | f 


“i & 
he yy = BL C, + dame, (My tQVP% iM.) +e + ae 


fn, 71 (2ite Mey 
ie (M,+2V"4, ) (Faia, + 28 


and then substitute the resulting inequality in 3.110. 
We note that the following terms appear on the right side 


of 3.110 after this substitution 





2 


pees Ley 3M M 
3.116 me LS + gman, (My +2VH,) + 3 +e 


+ 33 (My, +2\ iMG, ) nd, ta ESN 


i (4 +2V7,) ff a? ta Ys, 7 
es 
er Bs fay a %, ds —M Sf eds 


] 
ET tg V( d 
4 <A, qo) 


~ 
= 


where M3 and M are defined by the above inequality and 


as 


Beal7 Mm. = ae >on 
D(y, 


Since we are seeking an upper bound for the integral of the 
square of the conormal derivative over S, it is clear the 
unknown terms which appear in 3.116 may be neglected in the 
revised version of 3.110. In the original statement of the 
boundary value problem, the solution was assumed to be con- 
tinuously differentiable in VUS; if we make the stronger 
assumption that the continuous differentiability holds in 
VUSUD(y,); then it is easy to obtain an upper bound for 
3.116 which is less than zero. In what follows the stronger 
assumption concerning differentiability is made since, if 
this condition does not apply, then one may clearly 

neglect 3.116 in order to obtain the desired result. We 


define 
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3.118 B(xe) (7 Diy) 


and then write 


3.119 J fin Vd l= {Se ds +2SS t'vyds 








Dio) 
a+k) SS Pds eee ff o'ty ¢ 
i) D(4o) KG O ff a 34 J4 ds 
where 
3.120 M5 = max \ei } 
B(iyg) 
Next we set 
Bele) EF Meroe de = —_ 
mM. My 
and solve for A to obtain 
3.122 : ia + Vm? + See ) 
From the above it is clear that 
Boe 3 me Dom. Shay Yds —M it Yds 
| goad, Dex) ya ‘4 al. 
4 “Ge 
Kee 4,4, (-H,+\ Mi , Me Ma, | d t 
aM, : 


B (40) 





et 


The inequality on the right of 3.123 is strict since in 
Chapter II the auxiliary functions fi were chosen so that 


fin, had a positive minimum on S. 


it 
By combining the above, we are able to rewrite 3.110 


as follows 
3.124 BS )ds <4 Xm, SSF Pav 
+4 (My +2VH,)' Ye ng SSS F'dv +4 £ Lf 6,.a°%Gds 
4 Soin, (BB) ds + om, Sa Lam (oma JF Cs 
+A [Cras y' (34) AS + im UGA gammy th #2VF,) 
+i Ho + 2 (uy save)! (Ey, + BY} 


15, Sfp om Grade off emp onan 


" iff Fay aK, (y,-4)|Av = {Sim Ho > brs [2 ky (,-a)]ds f 


ath, (=H, -\2 Re ey f, W2de = 
Ms “ 
B( 





where a) is given by 3.115 and the computable constant B 


is such that 


3.125 gees Q 
when 
3.126 F=H=G=0. 


Inequality 3.12) is a bound for the unknown factor 
of the last term on the right of 3.8. Bounds for the two 
remaining terms on the right of that inequality may now be 
obtained using the procedures of the first section of Bina 
chapter. With ate known it is clear that 3.57 becomes 


Best 2q His. JS 124 VI?) - °F (F) 
ll 


12h Fat + a2n2h FY dy 


< We He np [4 IK (4o-4) lds t 6 


: | 


= C35 HF pup [2K (y,- “lds ~ ~ ne Hope k egal ds 
LG 26 op (any,) als = + 


The right side of 3.127 is computable and hence the desired 
bounds may be obtained in exactly the same manner as was 


employed in section 1 of this chapter. The results of these 





ee 


where Sy ae given by 3.115 and the computable constant B 


is such that 


A 
BE iiec B= 0 
when 
5.226 PeeweeaG = 0. 


Inequality 3.12), is a bound for the unknown factor 
of the last term on the right of 3.8. Bounds for the two 
remaining terms on the right of that inequality may now be 
obtained using the procedures of the first section of this 
chapter. With (8 known it is clear that 3.57 becomes 


3.127 Jom SSS (a G3 )- 02 5@) 
ape V Vv (4) 


t2 Y ah gy + 2n2e% V4} dy 


<S(F “He Dep WK (maids + 6 


5 
ee 
a IU S5 Hyp Eek Gy “lds ~ ~ Sn Fe Hoephka pds 
coe Co 2 a \ 
D(2) op ( Mo) a 
The right side of 3.127 is computable and hence the desired 


bounds may be obtained in exactly the same manner as was 


employed in section 1 of this chapter. The results of these 





19 


operations may be written as follows: 


_ Wes) ~ 2 GA 
<{SET@] T(t )tdv }* 1S [ser] (FF). 
phe, Gnajdve™ + 19, 
PIN FF ple dove) 


ae 


EMS F open g-pldvt {5 § 
: hata bY phe 


where the functions Nes and Th are given by 2.10 and 2.39, 


the constants K,, Ko, and K, by 3.110, jeecnand 35745 


equation 3.l27 defines \ and the constant o is given by 


a 
3.75. The computable bound ee obtained from 3.12) while 
X(p} is defined by 3.85. As in the problem of the previous 


section, we remark that if the function f is such that 
3.129 f(x,y,w,V w) = f(x,y) 


then the result lends itself to improvement by means of the 
Rayleigh-Ritz technique. 

In obtaining the conclusions of this section, the 
role of inequality 2.113 was basic. By referring to the 


derivation of that inequality it is clear that the symmetry 
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of the components of the matrix aij was essential; hence the 
consideration of a method for the symmetrization of an 
arbitrary problem is appropriate. The result is obtained by 


defining 
3.130 ats = S(ali 4 add) 
and then, since 

=) = @J 
3.132 a Togs a Ps 45 


almost everywhere for (x,y)EV, we may write 


= ‘4 
3.132 (ty i) aw, tary, 


rayy, elated) 4, 


saily ~aity + (ay) 


‘6 

where 3.13¢e is also valid almost everywhere in V. Due to 
the foregoing, the differential equation of the boundary 
value problem may be replaced almost everywhere by 

A} 


M 


3.133 (a** ree -255 5 = £(%, 4400) tay oe a 


F(x, 44, ‘ vv), 
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From 3.133 it is clear that f satisfies the same conditions 
as f and thus the original problem has been symmetrized as 
required. 

Before concluding our consideration of normally 
parabolic problems, we remark that when the differential 


operator has the form 

- as 2 
3.134 g(¥) (a Pega wy 257 ? 
the transformation given by 3.20 yields 


7 oe 
he - (4 + K)¥ -g : 


Sua3g 3(P) = (atl % yy 


J 
Hence, if the minimum of the function on over V is a strictly 
positive number, then the constants K; of either this section 
or section one are reduced by the amount of this minimun. 

In order to appreciate the improvement which results under 
such conditions, the reader is referred to the example 


presented in Chapter VI. 





CHAPTER IV 
DEGENERATE PROBLEMS 


1. The Dirichlet Problem at a Parabolic Point 

the problems which have been considered thus far have 
been normally parabolic in character. We now turn our 
attention to problems which are normally parabolic at certain 
points and degenerate at others. As the above heading 
implies, this section pertains to the derivation of a bound 


for the solution of a degenerate boundary value problem at 


a normally parabolic point. We seek to bound the solution 


of 
T(w)=(atJw 5), j-aw-edter (x,y) (x,y)EV 

H- W(x,0) = g(x) (x,y) € D(o) 
W(x,y) = R(x,y) (x,y)ES 


where w is continuously differentiable for (x,y) € VUSUD(y,) 
and is such that its second derivative in x and mixed 
derivative in x and y are continuous in the interior of a 
finite number of subregions, the sum of which is D(y) for 
every y satisfying o<xy<y,- In addition it is necessary 
to assume that 

lim w(x,y) = g(x) 


( sayy) 7 (X50) 
(xjy) © V 


ae 


Lim Ws (X,y) = g; (x) 1, = 1, 2 2 sae 
(x, y}—9(x, 0) 
Cxay) eV 
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The components of the symmetric matrix alJ satisfy the 
conditions imposed in the second section of the previous 
chapter. The functions q and § are piecewise continuously 


aifferentiable in y for (x,y)€V and g satisfies 
3 Z(x,y)20 (x,y) EV. 


Prior to stating the additional conditions which are to be 
satisfied by the functions Z and q, it is necessary to 
define a subregion E of V by 


Ll B= {(x,y)| B(x,y) = 0, (x,y) EV. 


We suppose that q has the form 


di; i 

° qg = E 2 

where the constants gf. e@i, 2s ss dL, satisfy 
Z be 

4.6 > é = / 
= 


and the functions qr iad, oP 2 fulfill certain 
conditions which will be prescribed in what follows. Let 
the coordinates xi, i=i1, 2... WN, be chosen to satisfy 
.7 min max S (xtext)? = R° 

(x, y) (x,y)éS ! 


where of y<y,. For N = 1 we suppose that 


1.8 a m = min}o3q> -a, (4R2)~* } 
E 
& 





with ECV defined by 


n.d B.={ (x53) | I(x,y)-pll<§ for some Dé€E with S>oh MV 


In }.9 the symbol || || denotes the Euclidean distance in 
the 2 (and later N + 1) dimensional space of the problem. 


If on the other hand N = 2, then q must satisfy 


Lig b ysmin{osryap =85 (URT) 4 L — aay 2 eo ¢ wh 

% 
where 

N 
= oO Ll 

ees in 
HO ry 4 (x x7) \s 
and 
gat R= 7 ae 


Finally for N 23 


1.8 Cc iiy-minjosrta; > a, (1-2)? } ile = L% Oo ° e dW 
& 


Since the techniques which are used to obtain solution 
bounds for N =1 or N =@ are very similar to those which 
are employed for N23, the method shall be presented in 
complete detail only for the latter case. Thus for N23 


we define m by 
i Ze 
4.12 O< “LE ih a, (N-2) =m¢l 


If N=1 or N = 2, similar constants are introduced. 
Previously we required that % be piecewise continu= 


ously differentiable in y for (x,y) € V3; in addition to this, 
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it is also necessary that a be continuous in V for every 
(x,y) € —~ where E ts the boundary of E, and 5% 70 on E N(V-Y). 

In .1 the boundary function f = f(x,y) 18 continuous 
in y for (x,y) EV and the first derivative of f with respest 
to y 1s continuous in the interior of a finite number of 
subintervals, the sum of which is the line (x,y) for every 
fixed x where (x,y)EV. In addition we suppose that 4 
is integrable and square integrable over V and that g and 
G», are integrable and square integrable over D(o). By 
definition it is clear that f is integrable and square 
integrable over D(o), D(y,)s S and V. 

It is now our aim to determine a bound for the 


solution of )|.1. As the heading of this section indicates, 


we assume that 


4.13 Z(x,y)>o (x,y) € D(y,) 


and that 2.9 is satisfied at the point p€ D(y,) at which a 
bound is desired. The arbitrary function ¥ = f (x,y) is 
chosen as in the previous chapter and the function Y (x,y) 


is introduced where 
4. 1h W(x,y) =w(x,y) = YP (x,y) 
We compute 


J(¥) = F(x,y) (x,y) & V 
h.1s Y (x,0) = G(x) (x,y) € D(o) 
i (x,y) * H(xyy) (xsy)ES8 
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Let 


5 oe. : 2 
1.16 J(¥) a (a as 3 qg ¥ * S¥ (Z VY), 
and then with ey given by 2.10 we note that 3.15 is valid 


for the present problem. From this it follows that 


4.17 IW Oxy "46) - X40) 


Clow fff v T(%>) dvi 
4732 V(4) P) 


+1 SS (24) 454% | Sf yt dst 


4-16 X(p)= P54.) — JS 32 d's 


Quen SSS FoF dv On, 26d5 - SSn,¥ EHS | 
47 fe VG 


Since F is bounded over V, it is clear that the volume 
integral ir: 1.18 exists; thus the desired bound is available 
when the terms on the right of ).17 are estimated. We shall 
first consider the integral of the square of the conormal 
derivative over § and then use this bound to estimate the 
unknown volume integral. For this purpose we make ths 
transformation of variables given by 3.20. The boundary 


value problem then bécomes 
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ae Ree =) oe YY a 
TP)= (a FD ~ % V-254 =F anplicy-g)]  Ggdev 
h.19 F(%0) = G typ (K4,) (x4) € D(2) 


Vosi)= 4 exp [k (y,-4)] («JES 


where 


4.20 qzuKZ+q., 


Due to the fact that ae is continuous in V for (x,y) é E 


and ..8 is satisfied, we may choose K>o so that 


~~ 2 — -| 
luteal, Cc A os Se ‘6 > My (4 -ne[ade") >-£4 (v-2)*| 
TP SL 
for N 2 3 and 
2 
. a: 1 O¢ =, 360 
wee min (53KZ = = ==) > = (l-m 
where 
a = AY = e I — Is 
b.23 q = Kg+ é“a, 2-é (Be otag) Ena 


In the above we have made use of the fact that 28 = 0 


on E. Similarly the constant K7o is chosen so that 

eet = mintq-5 2yot = 

ya fi = minyq-5BK >-a (4R™) 303 for N=1 
V 


and 
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4.21 b = “gre Oy act (4, -K2[ad ¢ = )>-a,(¢ey' 
T2120 and we2 


With K>o as chosen above we have from 2.113 


I 2h (wads < SS B5)As 
Me Be (6 +9 + te) Sif abd, 9, dv 
tg! SI Pd + im, SIF sep feng] dv 
= se (Ba 
where K, is an arbitrary positive number and 


a 
25 My 5 max | 8 > (24)* § 


1.26 My), = ” a} 


if 


In order to estimate the unknown terms on the right of h.2h 
we observe that for any set of N functions gi, Lo) VS eee 
- N, which are piecewise continuously differentiable in x 


for (x,y)EV, it follows that 


4-27 fain] ds = Shas Fay 4 allay. y dy. 


By setting 


1.28 gi = (xt-x1} 





Bg 


in 4.27, it is easily shown that 


4.29 


ma, 8 (hit) ee Ce 
Mf ¢ dv < WS itn) Peds + Fh Sa 4%, a, 


Substitution of this inequality into 4.2) ytelds 


. 30 


Lf (3)'ds 
< XM +2 ( Gt Ee te +e) [a'tg ba G, dv 
then Sf (x*-x ‘Yay, H? bup |R kina) W'S 


+ A SSS Fou Bk U4y-a)]dv + coe + Je (ey dv. 


To bound the unknown terms on the right of 4.30 we write 


h.31 


Sf 34 FP dy = Lf FF ls ~ Sf Gr on (ake,) ds 
oF 

~ SIE F 34 Mp goa) Y + Wk Pr ou [tyr a)]dv 

+ fin, We anp entra |S #- KIS 354 Rady 


av og alt 
55 4s -% LG Y d's ~# Sfm 4H orp [2k (4,- lal s 
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14 36 “ovp (2hy,) as + ESSO 53 =$ Jy - EE: a 


Due to the assumptions which were made concerning the 
aifferentiability of the solution of h.1, it is clear that 
3.93, 3.95, 3-97, and 3.105 are valid for the problem 
presently under investigation. Substitution of these 
results into a transposed form of ||.31 gives after making 


use of the fact that ¥ = y on Dy, ) 


I, «32 Lfe(¥ dv + cus te ds oS Ge) eds} 
Go 


<z(M,+2V%,) Dyp (2 M4) a at ds t z2, a Eds 
s 1, 


: 


+( za i ‘fo a ray) Sifath, dv + Dn 


where we assume that (Mo + 2VM 10) is a strictly positive 
number. If this is not the case, obvious modifications in 
what follows will yield the same result. Inl.32 it is 


clear that 
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Mn, = - sf my HF on RKC g) obs + i CF expo (Ruy) ds 

t tL st, -k F) up ek Ca,~~9)] dv 

EAS FC erplergyds -& Sn 5 Hong Pxcq-a) Jas 

ty fi GA % Dap (2hap) ds +97 (RM) x") mH ouplendy als 
? rSS a, ( BE) sep le « ¢y,-a)'ds 

+(Mt,) S i-miy( Se Hug [i laa) a s 

roy Sen, lam Come dl (BES ong fancynd] as 


and 


Wh as ome 188 | 


We now show that due to the form of q and hence q it ta 
possible to choose fy so that the terms in the braces on 
the left of ).32 may be bounded from below by a computable 
constant. On the surface D(y,) we have by means of the 


divergence theorem 
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4.35 4m y Ie 
- 


is 4, “wd A Y Wal s 


where Bly.) is the boundary of D(y,)- For N23 the form of 
1.35 leads us to define 

1 i 
le. 36 g cx xr) = 


and then for every I obtain from .35 


4.37 SS 2, a Ale 
D(a) 


fait YY. 


<< (x*ox" az’ ~ y 10 A 
(N-2)~] 3 af aT N-2)K - “I] DG.) 2! ot 


where 4 is an arbitrary positive number satisfying 


38 ” 
3 ae ae 


We choose & so that the coefficient of the last term on 
the right of }}.37 is minimized. The above equation thus 


becomes 
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Lwig8 ff Peis 
¢f ‘A\ oh ul 
re J 6 )A, m.¥ dL + ay (v-2)'1, 4 ffi ay de 
For later use we observe that it is also possible to show that 
4.39 ¢ ISS Pin Poly 
gh J a 4\ 2? “el alas y oy 
SNR (as nn y S + Wavy, fifa" a 


} 
For N = 1 we have 


aE RS 2 
38 @ S ytds 2 (x-x) 4 + 18 Say yy ds 
D(y,) Kx, ? D4) 
and 
n a > 2 
4.39 a (SPdv 4 £1 um ds + + SSal'y y dv 
0 x 9K 
Finally for N = 2 we may write 
: 8 b 4 Dui 6 f (ret yet a iz 4 th: At | 
4.3 am , d i) - em vad aye f 44 ds 


and 


4.39 b Se ng'dv € LL ot “Ay gm, 8 ods + BLTa' EG, 


0 
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We are now in a position to obtain a lower bound for the 
terms in the braces on the left of }.32. The details shall 
be worked out only for N23. Similar results easily follow 


in the other two cases. By means of the above and 1.21 ¢ 


oar 4 vd > LS eax Yands mee ds 
go 


‘ _ (aE) Sons af mH Al aim Sfalte 


fe D(A) 4 tia 


We notice that the line integral on the right of 1.0 is 
actually a volume integral over an (N-1) dimensional 
domain; and, since N23, it is clear that the integral 


exists. In .32 we choose 
2 


7 _N as 
Hed = (Lm) Fe >0 


where the strict inequality on the right follows from l.12. 


Substitution of h.41 into h.32 produces the inequality 


2 Lis a(S Jey <5 (M,+2VM, ) amp (an) SS (5) ds 


+2R eye Q sim) Se e(fe + 28 18 2K wig) a & a 
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LSS 224 Ge. i a A 
Tt ee | (ora), Sy ‘mn Hd 2 


t (ina) Bar J (x-x") my, Holl + hy 
~ BH) 


In deriving the above, an inequality analogous to ).29 over 


D(y,)} is used to bound the integral of the square of YW over 


D(y,)- 
So that ).30 may yield a non-trivial bound for the 


integral of the square of the conormal derivative, we choose 


OC, as follows: 
hh 3 O, = 5m5(Mo+2 VM) “*exp (-2Ky,) 
Inequality ||.30 then gives 
, ot)? Ps 2M, 2M 
4 Ly Sf (33) ds < orig 1G + Go + G3 (My « ava) orp (are) 


2 2 2 
4 Py RY (<e y 2k Ms ~/ 
BME so (ft REO) e203, 


where 





Ws 
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= 2H + mea Slt), Paep len cgy- ad]ds 


t i i , J F aup i d]dv ig 2(My+2VAy,) 2p (-21. 42 Iv’, (a) 


SfFads th - Soa, w= )f G's iginy Wo 


+(l-m) ee ‘. Gx’) mH ALS 


The unknown term appearing on the right of h.lh is estimated 


by using the divergence theorem tc write 


146 


Sf al" ¢ ha ty av = — SIS PtP) dy 


a oY ene mee 
NY Ds ds - ING dv - 351 2 Ps 


-E£ (fn 2F'ds +4[l2ePols rL (Se Vy 
ca D(o) V “3 


: +o Sf ate ¢ dvs SS (x'n'Yn H'aep[2k(4,-a)ld s 
ta Sf Feepkingpldy + IS 3) ds 

3 H dep 4 K(4o-¢) | ds 
= Sha (ger rT) SS ag? mH Dep l2K Cape sas 
+m (150° . y dv i ze 6 ep ak) ds 


dL 1 





where C4 


me 


~RSSn,2 Hop lRK (ye ads +£(i- a) fate va al 


+ (im) Fe SS (x) Heap [ax(g-a)]ds 


and t are arbitrary positive numbers which shall 


be suitably chosen in what follows. To obtain the above 


it is clear that 4.21 c, ).22, ).29 and ).38 c have been 


used as necessary. Inequality l|.46 is now transposed so 


that the unknown term on the right of l.4h may be bounded 


47 


3(I-m ) Sff a’? a ¥ % dV 
cz SS SCS) ds ex LSS Hou [Hk (4-d)]d 
+( “ ‘SS ea')n teepfaneysd)]d 


+ (1 \*< J JJ Foe aK Cg ~4)] AV 





ch 
“26 ay N- Coico “x. Jag, Hovplak(g-)| ds 
+4 [fee'vep (any) ds — % [ng 2 Hauplin gry) s 
+ (ina) aR JJ (x'- x") m Heap 2Ku-¢)] 4s 


> 4 (3)ds +E SSH ceplt rigs d]as thy 
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where Mn is defined by the above equation and it is clear 


that we have set 


es (iam 


Since m1 due to 11.12, inequality 4.47 provides a non- 
trivial bound for the Integral of wD Y 5. By comparing 
Lek and h.l7 we are led to choose 


= (def 4, “1 


~ 
- 


then .li) becomes 


l, 50 IJ (4) <2, t In, (5) [Ie oypten ras 
(n,m = B 


where B ia computable and is auch that 


4.51 B = () 
if 
4.52 Pe HeG #0 


Tne above is a bound for the unknown factor of the 
second term on the right of 4.17. We next bound the volume 
integral on the right of that inequality in terma of 
computable constants and l.50. Let 


\.53 M(H) = (ald Wiad, aah + &y (BY ) 





Ne) 


and consider the function 4 given by 2.7). Since the lower 
order terms of J* and J are identical, we may deduce from 


Theorem 2.3 that the arbitrary constants of 2.7) may be chosen 


so that 

h.Sh -~-BomM and IH) <o for (x,y)EV 
and 

4.55 I°(F)<o 


for every (x,y)€ V such.,that o<£&< ¥o<7 $7: From the 
proof of Theorem 2.1 and the definition of Y and g, it 


follows that 


1, 56 Lorn. I H2y“ds=o 


By means of the divergence theorem and the differentiability 


of p we have 


4.37 Sim. SIS I(Y)-V TY dv -LLg vv 
47S td) ’ 
4 *Go 
Si fF BF Si dds — ine vd 


hi |Fz vd s 
Dio) 
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and by expansion of the last equation the following inequality 


is established 


l, 58 bom See yay) YH) dv <-2SIFa'ty.y, dy 
74, 
4d 


SIF - Si-m, 2 Hds + i Gds+b 


cam + SS {pe as 2% f Wd I Fee ds+B- Y 


where J is obtained from L.u7 after application of ).50 
and M is given by .54. By an argument similar to that 


which was used to establish 3.64, we are able to write 
1.59 ASF eT H dy - SS eT) dv <\) 


and finally from equations analogous to 3.66, 3.67, 3.68 


and 3.59 we have 


l..60 Pion. US SSS pT Cm) dv 


<t SSS bra] (8) af y(t { : Serre Fe y 3% 


+9, + UES me Fv}. 


The existence of the integrals on the right of .60 is 
established by a discussion similar to that of the previous 


chapter. 





LO 


Inequalities ).51 and .60 give computable bounds 
for the unknown terms on the right of 1.17 and hence the 


desired pointwise bound is obtained from the inequality 


4.62 IW (x, 4.) - X(p)\< {Sere (ya 
Serer] Fry +0 State] Fray > 
[Sutrdst {Bt 


where the functions 3 and Gh are given by 2.10 and 2.7) 
and the computable constants B and 9 are obtained from 
1.50 and 4.58. As in the previous problems, we observe 
that the right side of ).61 is identically equal to zero if 
F=G=He=0. 

During the investigation of the problem of this 
section, an attempt was made to generalize the method by 
omitting the requirement that 23 be continuous in V for 
(x,y) € E. Such a generalization has not as yet been achieved; 
however, if the domain V is restricted as indicated below, 
then the desired bound is available. 

We suppose that the given boundary value problem is 
identical to that stated at the beginning of this except 
that = need not be continuous in V for (x,y)€E. In 
addition, however, we suppose that V is a cylinder (that 


is; D(o) = D(y) for o Sy<y,)- It is evident that equations 





LOZ 
4.14 to 4.20 are valid. We choose K>o so that 


nw on a 2 = 
1.62 fiysmin{os rtd, > ~fa, (2) } I =e fe e/ 


for N23 where it is clear that K>o may also be similarly 
defined for N =1 and N =e. In what follows we shall 
consider in detail only the case N2 3 since similar 
techniques would give the same results in the other two 
cases. 

The integral of the square of the conormal derivative 
of ¥ is obtained by botinding the unknown expressions on 
the right of ..30 in terms of the given boundary data. From 
the divergence theorem and the differentiability of W and 
the coefficients of the differential operator, we have for 


an arbitrary constant "a" 


t~ 


1.63 ier S(%) sup faly,-4)] dv 
o WF 2 AtZ2k 
bay a — 9 earl 


(an) Sf YF op (arr) (4-4) AV 


—~ 2 
SF By plan ised] dv < 





(E+ St) Ma ee, oop htardpv -4 Mf 0b, Ga 
Hl a%G, 6, oup[laran)gds + $10 (3) ds 
+ 2E. HE O44 cle exp [2 (at2k) Cy, ~a)] ds = ri ge v ‘ds 
raf J 4 ee Mp) laren) |ds = a9 y Srp [aty.-a) | al 
red 54 Pipher adv - 5 iie(35)) mp b-a)| AV. 


The above inequality is transposed and use is made of .29, 


4.39 c, 4.62 and Schwarz's inequality to obtain 
— {2 fe fam - She wa bra" HY 
Ma-ES Sah yu + ie ag 
cof {S(t +i If] 35 -ki| deple (ate) yr8)]d' 
rey a6, G,, wyp|(ar2r)s,]ds + if 4 G inp | aren) ds 
ig sf Eds 4 if er pl (ark), ds 
+E SST RE - (24) exp eCard (ye gl] Av. 





Lou 


In .30 and the last inequality, we set 











m 
h.65 ¢, = = 
4.66 c. = 

N@a, _ 
ie 7 c, = rf (1 - m) 
and 

Mz 2M, Re ee 2 a a 

l,.68 a= (1am) [ie + Gd + hie, tam (Gt Fe tae +E 


Substitution of the above into l|.30 yields 


une (Sods <2H + mht JS 2 I, Hin x (4,~4)] A 
mm, SSS Poaplegea)|dy + ISLS ki oul @renty,-d)] ds 
+ Sfa%e, 6, omplareny dds + 1G Gap arena.) As 

rma, (=) ff Fas rdf 6Foxp[(arera.] ds 

USE Cmdr] op rang] dv = B 


where a is given by )}.68 and K>o is chosen so that l.62 

is satisfied. Since 3 is a computable constant, the 
desired bound for the integral of the square of the conormal 
derivative has been obtained. A bound for the unknown 


volume integral in l.17 is again given by .60, and thus the 
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desired pointwise bound is available. 


2. The Dirichlet Problem at a Degenerate Point 


We seek to obtain a bound for the solution of ).1 at 


a degenerate point; that is, at a point p which is such that 


4.70 Z(p) = 0 


At a degenerate point the author was unable to determine a 
parametrix with properties analogous to those of ae at a 
normally parabolic point. Because of this reason, additional 
conditions had to be placed on the boundary data, the shape 
of the domain and the region of validity of the differential 
operator. 

Let the problem for which a bound for the solution is 
desired be that which is proposed at the beginning of the 
last section (specifically, from the heading to the paragraph 
immediately preceding equation }.13) with ).70 replacing .13. 
In addition we suppose that V is a cylindrical domain, L 
is twice piecewise continuously differentiable in y, and the 
differential equation is satisfied on D(o). Moreover we 
assume that the derivative with respect to y of the conormal 
derivative of w is continuous in the interior of a finite 


number of regions the sum of which is S and that the fol- 


lowing are satisfied 


eae en _ fh 
p> PEB(o)? ) ba (A) 


P, E Dio) Pre S 
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inet 2 Z(x,0)>0o 
h.73 Lim 2275) ee o) 
ay vy 
D7 ( x50) 
peév 
and 
4.7 2(X,y,) 


for some open set D'cd(y.) for which p is an interior point. 
To shorten the aoe we consider only the case 

N23 since for N = 1 or N = 2 similar results are easily 

obtained. As in the previous sections, we obtain the 

desired bound by introducing a function Y = Y (x,y) with 

the same differentiability as w and then with Y given by 

il we compute the system .15. Next we define the function 

y given by 3.20 and compute l|.19. The constant Ko is 


chosen so that 4.21 ¢ ana 


2] 
lig, 22 aa 
h.75 (SKZ + 5 327 “(ae 


are satisfied. The desired bound is now obtained by using 


(6.8) of /21 | to write 


L.76 ¥(p) = pou PAY) - Yar ») 5 ds 


giv sy - 2 3} de 
+o ~ 





uO] 


where 


4-77 (L(Y) E (a"*¥,,),, 4 Y 


os ty \ 
(a an me aa 
and a is given by (l.2) of [21]. The unknown terms on 
the right of l|.75 are the first, second and fourth. For 
the fourth term consider that 
2 a, te sy 4, 
1.78 $f, sida<] sa duh K: (Soy als ¢* 
Hd; ) BCq) BC) 
Using Hormander's Identity and the same techniques as were 


used to obtain 2.113, it is clear that 


4.79 $ (dls H+ % Se*y atv). ds 
Bé4,) : D(a) f w$ 
RC, | 
where 
— “yb —| i A 
80 H z wn. g feng) (4, te at ) AL 
+ 7, nt fem, stp 
B (do) 
mas ’ -1_k 
1.81 m, = By.) fon, 


ana fl = fi (x) are any set of N piecewise continuously 
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differentiable functions which are such that fay is 


bounded and has a positive minimum on B(y.)}- From h.32 
after application of .50 and ).h7 or from .6h with the 
unknown term on the right bounded by }.69, it is evident 


that one may obtain 


—Sae 


82 SS atyedss DB 
D(4.) A ~¢ 


where Lf is a computable constant. Inequality ).82 bounds 


the last term on the right of ).79; for the second term we 


write 
1.83 2s ¢* (aby, Je as 
D(4.) 
all gt a Wds + fle eh tsp 3 ds 
Dy) 
f 
+o FR y, ds 
<(Ge 7 ) fsa" ia "4 4. AS r My ek 
2 
+ {fz (38) ds + SJ F*ds 
Did) (qu) 
where 
_ 2 
a6 Mi, = max 4 >( ty 
— "1 met a 
M 
4.85 i, = 7 a Ley” } 


and 
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L,. 86 = max (a Bi} 


D(Yo) 


By using an inequality over D(y,) which is analogous to ).29, 


we may rewrite l.:79 as 
87 § (3440 < 
5 (ay) 
(C+ +F a, yak) By mat $e xn dl 
t& SS eds 3 2(3 
mi = Dy) + > Diy (Nv 


Z a avy 
mn, tn, {Se Sq) dv 


For the last term on the right we use the divergence theorem 


as follows: 


ait 


2 ip Bea ey 3B 2 gy 


a 4 


ier Shs 5 


V 


S34 By Bs - M7 BEY Voki 2 (Sq )ds 


resp2($ st yds - x SSS OS 5) aya 
DY 


Since 





wee PS ee ee ee eee 


LO 


~ 


‘4 
Sa Oy c Me ate bt ty 
ne vy, By dy Seg BY, + Fe mi Oy 


From .4.7 after application of 4.50 or from l.6l with the 


unknown term bounded by .69, one has 
4.91 Sate ¥, dv = poe 
V 


By transposing |.88, and integrating by parts, we obtain 


he sf a( pares Egy MR Hl Bee alt Sd 


Pres (2) dO 4 xh, af [rm sg] de 


- 
ae. # 


+H Ve -x'\n Faw + Be) opp akKe.-#)] 
+, JSS KE 4 35] oxplen(y-)] Ay - 


where 


eae 


6 


4.93 Dy = >a — {-kH4 re exo (21, ) Al 
U 
te Sli 2K 4 Sait aph ty.-a) ds +2 


3 ei, (Ea) SV a agin fates Be Vp ence 


i) 


ELLE LNG) ge FE empty.) de 


D(o) 


+(I- m) oP JEG xan, [-K#+ Sh Uap knty- 4) ds, 


In .92 we choose the strictly positive numbers o., and ec 
so that the coefficient of the second term on the left is 


zero. By comparing ||.87 and 1.92 we are led to define 
Ly a =i 
: = ene 
is 


By fixing the arbitrary positive constants as indicated 


above, it is clear that one easily obtains 


>¥ 4 Vy 
h.95 € (Sy dl < RP 
B40) 
where 8 is a computable constant. Inequalities 4.95 and 


41.78 provide the desired estimate for the last term on the 
right of |.75. The second term on the right of that 
equation is bounded by using exactly the same procedures as 
Payne and Weinberger employed in [21]. Let (p} be defined 
by UO aor [21] ; it is then clear that le and & 





Lie 


X 


may be chosen so that 


4.96 ae for (x,y) € D(y,) 
and 

297 a* F) <0 for (x,y) € D(y,) 
where 

4-98 CK) = (aE) 2 BF 


From the divergence theorem and the definition of 
4.99 if 
Sf {ha (%*)-¥ a) ds Ley’ 12a ty, yds 
Di, DG) 
=- {f va @ds+ 2 vFaods 
D(40) Nao) 


2 mee yD. 2 SP — @ 
ci Hdl + B 3 sy 4 L+ 2h d= 


and thus 


Se eer 5 
100 uff. Ya (F) ds} 
s{ Sfear@y'F* gcytd st” 
Dio) 


ae “Ig 2 | SA, 
+Y {Sa Ne *ayrds +09" 
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To bound the unknown term in the above expression we make 


use of 1.82 as follows 


Pot (co 
4.101 Sf wads « & P(x x Jn, H'd's + vo 
Dé4o) B40) 











and then obtain 


102 fede l'F*aw* ds 
D(-4y) 


OPy\2 - -| 
SRM, {f 2 (Sy/4s t 2 bal) FFs 


<4 MH My Oy +4 M4, Mn, +2 fercay'# rds 


where 

4% See 
4 «103 M, 6 nme, Lf (hr Gr g | 
10h = max {,,8K?} 


"17 D(y,)- 7 
and from 1.92 after application of 1.95 
Sie ye 
4.105 Sfz (Si )ds ch 
Ya, 6 


By means of Schwarzts inequality 


I, 106 Ny ¥ ACID) as 


h = : / 
<{if- vate) {Seta y acreyas 
Dio) Dido) 
In the above the first factor on the right is bounded by 


4.100 and .102 while the second may be computed due to 





114 


Gia b) of [21] . To estimate the first term on the right of 
4.75 we consider 


I 107 SJ raty) ds 
Dlda) 
. SY 
i etpsqds + If ry F dv 


SRM ny + ® Mg Ine + HG nF AV 
where 


4.108 = max ere 
Mg D(y,)-D" Pp } 


Transposing 4.75 and taking the absolute value yields 


l, 109 IW(p) - ¥(p)- ir pfds -9 H Sy aye | 
B(Ap) 


$2 M5 Ma My, +2 Mg Mg +f vac) els | 


he ves at 2 


where 11.106 and 1.78 bound the last two expressions in terms 
of computable constants. 

It is the belief of the author that the conditions 
imposed on the boundary value problem of this section are 
excessive. In future investigations, efforts shall be 


made to produce pointwise bounds at a degenerate point 





1i5 


without imposing the rather restrictive conditions which 
were necessary in order that the above techniques yield the 


desired result. 


3. The Mixed Problem at a Parabolic Point 
Thus far in this chapter we have always required that 
the matrix aid be positive definite and that % be @a non- 


negative function. Now we suppose that % is strictly positive 


but reduce the conditions on the matrix by requiring that a*d 


be positive semi-definite. 
Let V be an (N + 1) dimensional domain and let the 
notation of the problem be as indicated in equations 2.1 to 


2.4. In place of condition 2.5 we now require that 


Lee © min een 
The 
We shall derive pointwise bounds for the solution of 
= (y Kt ey 
J(w)= ( Wy? Oy ~ tOsqw) anne 


tepals) W(x,0) =4 0%) (x,y) € D(o) 
(x,yJES 
Ey = L x, 9) 
where clearly it is necessary that 


fate £ (x,y) =o 


for (x,y) ¢ S, with 
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eee 3 a | (x,y) | (x,y) ES, Yn? = 0} 


The solution function w, the coefficients of the differential 
operator, and the boundary data have the same differenti- 
ability and integrability as the analogous functions of the 
mixed problem of Chapter III. As in the previous problem, the 


function ff satisfies a Lipschitz condition in w, that is, 


Iiveralialiee If (x,y,w,)-£(x,y,W) |<, |w, =a 
As the title implies, the solution is to be bounded at a 
parabolic point péD(y,}. Thus atJ is positive definite so 
that a, 3 (Pp) exists. In addition to the above it is 
necessary to assume that conditions 2.8, 2.9, and 3.5 are 
satisfied at p. 

We use the notation and techniques of the previous 
mixed problem to obtain the following equation which is 


analogous to 3.16 


4115 LWoey4o- XC) S| Lim SSS ¥ F(8,) dv 
go) X Cp 440 Vy) ( ») | 
we SIS % LF saw) - Foo g, 0). aVI 

HL SurdsT HLS (Bp +, 2 xp)Mds 3” 


where 





ay 


4.116 A (p) = (x,, 44) + {Jt Yo d's 


+ Speeds - Pee eahye 


Go 7Ge v4) 
The right side of 1.116 is computable by the arguments 
associated with equations 3.18 and 3.193; and hence the 
desired bound is available when the unknown expressions on 
the right of 1.115 have been estimated. To estimate the 
boundary integral the change of variables given by 3.20 to 
3.22 is made and then by means of the divergence theorem it 


> 
= 


is clear that 


4.117 IJ Pids <2 (m my tz JJ ods 
+ ay i) Loup [2 il-\jel s+ Ck, +A, +4 +2 My) LI Pdv 
Ee LJ F* eryplan (s.-4)]dv G zife 6 ‘up (244,) ds ¢ 


where 


7m.>0o 


Liggan.8 min 3 = m,?m, 


S 
In order to obtain a non-trivial bound for the left side 


of 1.117 we set 
4.119 7 = $(mzmg) 


and 





Tete) 
=-{ Aw ee 
4-120 Kenny! (4%, +& +4M,)24, 


Due to the above we may write 
4.121 SS Yds <4 (9, om) SS Noyp 2h, 45-4) ds 
th Canzone)! SIS Foyle, (y,-a)] dv 
RC, acy! Sf ee Mapa (Ri, q,) AS 


= 6 


In order to bound the volume integrals of 1.115 we 
introduce the function 7 given by 2.39. Even though the 
differential operator of this section does not satisfy the 
conditions of Theorem 2.2, an inspection of the proof 
reveals that the conclusions of the theorem are valid for 
the present problem. First we observe that constants 
Po? SS > and OY associated with equations 2.36 and 
2.37 exist since J is normally parabolic at the point p 
and thus alJ is positive definite in a small neighborhood 
of p due to 2.8. Because % is strictly positive in V 
it easily follows by remarks similar to those which were 


employed in the proof of Theorem 2.2 that the constants 





119 


aw and # of 2.39 may be chosen so that 


fhetee Dear o, G<o (x,y) EV 
and 
4.123 ms 6, 7 <0, (x,y)E{xsy | (xs y)EV, 0< € 4 yo-y}. 


Due to the foregoing, equations 3.9 to 3.75 are 


valid for the present problem with 
a ek 
Het M5 = M =. 


The bounds for the unknown volume integrals of ).115 are 
hence available by means of the appropriate equations of 
Chapter III. After apvlication of these results one 


obtains 


6125 IWor ya) ~ 2 ~)| 
<{ SIE Fe dv}? 
CUES I (F Ff on fer g,-p] dvi 
19, + IESE (Fel cap feretye- ad] dv > 
HIS CF epplanca-aldvd {92 


HSS OS + my emp fds} oe 
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where op is given by 2.10 and ae defined by 2.39 subject 
to the remarks of this section. Using the results of 

Chapter III we obtain that Ky and K3 are given by 3.62 and 
3.7ha while the computable constants , and 9 are defined 
by 3.57 and 3.75a. The bound ® is determined by the methods 
of this section and is given explicitly by 4.121. 





CHAPTER V 
THE ELLIPTIC PROBLEM 


l. The Neumann Problem 

In [21], Payne and Weinberger obtained pointwise 
bounds for the solution of second order interior elliptic 
differential equations with Dirichlet or mixed boundary 
data. The occurrence of the zero eigenvalue of the free 
membrane problem nRomented the authors from extending their 
results to Neumann boundary value problems. By following 
closely the procedures which were used to obtain a bound 
for the solution of a normally parabolic mixed problems,. 
we obtain such bounds for an elliptic problem of the same 
general type. 

In what follows we shall use the notation and many 


of the results of section 6 of [21] . We seek pointwise 


bounds for the solution of 


B(w) = (ally ra) jew -qw = f(x) x€D 
9 3 


i yi 
2. Q(x) x€B 


where D is a domain of Euclidean N-space with boundary B. 
In order that the divergence theorem be applicable it is 


necessary to assume that w is continuous in D, continuously 


re 
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differentiable in D and is such that its second derivative 
is continuous in the interior of a finite number of 
subregion, the sum of which is D. In addition it is 
necessary to assume that the solution of the boundary value 
problem is such that 
5.2 lim atl (p)w, (p)- = atJ(x)w,.(x) {= 1, 2x see 
p77 EB J J 
pre 
We suppose that the components of the matrix aid and the Lt 
are piecewise continuously differentiable in D while q(x) 
is bounded and strictly positive in the same domain. As in 


(6.3) of [21] it is necessary that 


NV ’ 
5.3 oc LID Brqmls alte s, + b'¢ @, tea, 


<h CS BE + 9) s*| 


for all real numbers a ae ce) where equality holds on 
the left if and only if all the numbers are zero. Let 
5.4 min q(x) =m, >0 

D 


where my, is positive since q(x) is strictly positive by 
assumption. From (6.8) of [21] we have 


5.5 m6) - (p) = {P41 B(w-4) - (w-e) Bp) Ay 


baal 


+§ Cue) (Sof rm, f,)- rn, 5 (w-) $d 
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where 
5.6 B( lp} = (atip ) CBE) wap 
DoJ et 


with ES given by (l.2) of [21] and W defined as in the 
second paragraph of section h of that reference. We apply 
Schwarz's inequality to the first boundary term in 5.5 and 


6.7 $y (Srht by, Tbs {Gwe} *] 6B erat hols}? 


where, as in the previous parts of this thesis, we have set 


5.8 Y =w-uy, 


Let tt (x) eats pN (x) be the set of functions introduced in 
sections 2 and 3 of the work by Payne and Weinberger, that 
is, the set tt (x) is piecewise continuously differentiable 
in D and is such that fone is bounded and has a positive 


minimum on B. For later use, we set 


5.9 m5 = min t*n 


B k 


and then bound the unknown factor on the right of 5.7 by 


writing 


5.10 Suds 5 $ f'n ds = Lr # Dd v 


=m, SSH PP dy + am, SYP E: dV 


on 


< Fm, CPHL EY + aa ; 2 Swed 
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where 
eo M, = ma fs { 
5 1 a 4 
_M 
5.12 iM = max{> £2 | 
2 D / i 


and ee is an arbitrary positive number. By means of the 
boundary value problem and the divergence theorem, it is 


evident that 


kan.3 Svedv=-Sfya'ty, du +$ 42 ds 
D . Dann C 


Seb adv Sah cl 


where 
5.1 F(x) = f(x)-B( ) 
and 

5.15 L(x) = £(x) - 


we 


By means of 5.3 and a transposed form of 5.13, the 
following holds for an arbitrary positive numbers € 5 


and f 3 
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5.26 I ISE WE + (om $e) Phdv 
<ffaiy Ady, dv + LO %, vale Sig g)y*dv 


<b Sf ew +39 Bal ‘ a $ Lo ds 
: dD ts e 8 


In order to obtain a non-trivial bound we choose 
ol 2 —- m + 
5.17 ; ae 


and hence 


5-18 (f> y+ y| al V 


Sim SL Fe) “dy + E feds + sea § $ Lo) “ds. 


So that the last equation may be employed to bound 5.10 


we rewrite the previous equation as follows 
5.19 Std < <i Gam) IP vid + tlie I> ud dv 


and then choose C. so that the coefficients on the right 


are equal. We solve the equation 


e + OM €. = mM. =O 


Seo ee 1° Mo 


and obtain 
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Bael cS = ie |-21, +\ he + Bm, M | 


Since é, is to be a strictly positive number, the positive 


square root is chosen and thus 


—_ V2 mate 


To shorten the notation we let 


Eee C4 


meee |e 2 
5.23 M3 = m, My 3 ( ML, + M, + em, M.) + | 


and bound 5.19 with 5.18 as follows 


e _ ain) 2 M3 2 ™M 2 
5.2h ( 25, fuids <pepe Sf Pdy ‘zed, 2 Uds | 
Since a non-trival bound of the left side of the last 


equation is desired we choose the arbitrary number c. 


in an appropriate manner: 


ko 
5.25 C, “i 


Substitution of 5.25 into 5.2h yields 
2 
5.26 § 4 <n ts Sf Fay 4 i gutds = © 


which gives an estimate for the unknown boundary term of 
5.5 when this term is first bounded by 5.7. 

To bound the unknown volume term of 5.5 the following 
cases are considered. For N = 2 or N = 3 the results of 
pages 562 and 563 of [21] are applicable since from 5.18 


and 5.26 we have 





Lah 


5.27 x y* dy < Ema Sf Pedut & [BPG ds}*= 8 


A 


OT” 
For N >3 a function 7? is defined by (.18) of [21] so that 


28 “Poy x€D 
and 
5.29 B(@ ) <0 Gen, 


From the divergence theorem, it can be shown that 


5.30 S19 v B( Y)- Y BF )tdv< 


< 2mh O& + ERI Sets} 





+4 B = y 
where 
5.31 ", = max { t + b'n,® | 
Hence 
D+ 32 {SSP a@)dv b* 
D 


< i SSeBery FF dvs 


{SER FF dy + QFE. 
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The desired result is now obtained by applying Schwarz's 
inequality to the unknown volume integral which we are 


seeking to bound. We obtain 
5.33 Sa on avi Sacwww*dvP {I EEay [scqy)av} 


and note that the first factor on the right is bounded by.5. 32 
while the second is computable due to (las) of eal. 


From the above it is clear that 


5+ 3h IW(p)-2.(p) |. 


< e}"49 (bn Yds} 7 | S{yB (lp) dv | 


where 6 is given by 5.26, the second term on the right is 


bounded by 5.33 and 


5.35 “(p= - Ppl ds + SrpFdv + Wp) . 


2. A Generalized Dirichlet Problem 

In this section a pointwise bound for the solution 
of a Dirichlet problem shall be obtained where the given 
differential operator is such that the results of [21] are 
not applicable. 

Since the problem and the techniques which are to be 
employed are very similar to those which were used for the 


Degenerate Dirichlet Problem, the notation of [21] shall be 
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altered so that it more closely follows that which was 
previously developed in this paper. 


We seek &@ pointwise bound for the solution of 


ieee) = (ally -) 2s-qw = f(x) xED 
gts] 


5.36 

w(x} = h(x} x €B 
where D and B are as defined in the first Section of this 
chapter. The solution function w is assumed to be continu- 
ously differentiable in D and to have a continuous second 
derivative in the interior of a finite number of regions the 
sum of which is D. The symmetric matrix aid is piecewise 
continuously differentiable and satisfies 

Vv N 

5437 o<a2 Saves <a> eh 
for all real numbers (& 4, -.- - S ay where equality holds 
on the left if and only if all the numbers are zero. The 
bounded integrable function q = q(x} is very similar to the 
analogous function of the Degenerate Dirichlet Problem; that 


is, q has the form 


LD 
5.38 = Ely, = 2 74 
T=} r 
where 
Ld 
5.39 2. gta] 


nl 


a 


and for N23 
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‘ a= iy i 2 
50 mn, myn. {0 5 Ae gr? - 79 (W-2) } T=/,2,...8 
Conditions similar to .8 a and ).8 b may easily be obtained 
for N= 1 and N=e. To shorten the argument, we consider 
in complete detail only the case N23. As for the degenerate 


problem, we define m as 


51 0< -h eta [a (n-2)3| <1. 
So that the divergence theorem and Schwarz's inequality are 
applicable we assume that f is integrable and square integrable 
over D. | 

To obtain the desired bound for the solution at an 
interior point p€D we choose an arbitrary function / = Y(x) 
which is such that J(UY) and YW (x) approximate f(x) and h(x) 


as appropriate. We further require that Y be twice piecewise 


continuously differentiable in D. We introduce the function 


5.42 ¥ (x) = w(x) - Y (x) 
and compute 

J(¥) = F(x) x €D 
~ an): = Hie) x€B, 


From (6.8) of [21] we have 


5 hh W(p) - L(p) 


“SG T(Y-¥ TCL) bd $ sls I, Sohds 





rou 


where 


ait 


aN? TOPE) = (TP ads 


p 


with rp given by (1.2) of the work by Payne and Weinberger. 
The unknown terms on the right of 5.lh are the second and 
fourth. The final term is the first to be estimated. With 
B replacing B(y,) and D replacing D(y,) we note that 
equations |}.78 to l}.81 are valid for the present problem. 


Next we use 





5446 LS y"dv < ‘Ho ven Wds +E SSaty, w, dy 
to write 
2°47 - are S H + & (¢, Ae + va, fa" , bi av 


where 
5.8 m3 cs - ne fn, 
5.49 Ro = min max y (xt.x1)¢ 
xéD xB | 
ad 
5.50 Mi, = max 2 (f+) 
- D j 


5.51 My, = max a2} 
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and, 4s indicated above, K is given by .80. To determine 
a bound for the unknown term on the right of 5.7, the 


divergence theorem is used as follows 


5.52 Yes dv =~ SS yaw dy + $4 at ds Sg vay 


df 
es As oY a ae R i Av 2% 
x SS, a Ys dv + PH sv ds at (Ah, rE) 8G “Ke m HAs 


where we have use of an appropriate form of }.40. Trans- 


position of 5.52 yields 


ac, R° 


5.53 (1- m- wa, ) SJ a" 4, ¥, alv 


Q 
za iyo 
< ade oh, (2) $e) arm, Hd tte i) F’ aly 


XX SYA de a aL ‘og 2 
+2 § Grd t 2M, yy ds + 7 Gx dm, HH ALG 
for any positive numbers eo and L.,- Since we are seeking 
a bound ror the integral of the square of the conormal 
derivative over B we are led to choose the arbitrary 
constants in the last inequality as follows 


tae 
5.5 A, = (1-m) [Re 





= 


and 
ae =: 
je eee, ae 
5.55 dy = pm) (6, + > a | 


Substitution of these results into 5.l,7 yields 


ay A 
5.56 F(S)ds < B 
B 
where B is an appropriate computable constant which is such 
that 
eT Beo - 
ae 
5.58 F=H=0 


For N = 2 or N = 3 we use the results of pages 562 
and 563 of [21] to estimate the unknown volume term in 5.hl. 
These results are directly applicable since 5.53 provides a 


bound for 


N 
5.59 > v* dv. 
pa A 
If N>3 then # is defined by (4.18) of [21] so that 


5.60 $<! J*(F) <0 x€éD 


2S 


where 


Ceol J(F) = (allF ) -2qF | 
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The divergerice theorem may be used to establish that 


5.62 Sey TY)- Yet dy 


CAM dy +2 {pt { §F*u* als} 


and, since the right side of this equation is computable, 
it is evident that the deste’ pointwise bound is available 
using the techniques and the results of the paper by Payne 
and Weinberger. In the above equation B is given by 5.56 
while 4 is the bound for the integral of ald Ys Y 
which is obtained from 5.53 after application of 5.56. 





CHAPTER VI 
EXAMPLE 


Ll. Explicit Non-Linear Mixed Problem 

In this section a bound shall be obtained for the 
solution of a non-linear boundary value problem which is of 
the type considered in the first section of Chapter III. 

We suppose that V is a three dimensional domain 
bounded by the hyperplanes y =o and y = 1 and by the surface 


whose equation is 
6.1 ele + ao = ]. 


The boundary value problem is explicitly stated as follows: 


mw, Uw Ww =e 
DM) SE Se (at lw}) -42°+2 
6.2 K& yyEV - P< WV < 00 


W(jo)=0 Gy HVEDO) | We) AES 
where 
6.3 A >O. 


From the above we obtain 


Te 
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a 
eae vm, =! 
a 
2 mi =i 
6.) 
uh. 
a a, = V 
M, =o ae / 


By inspection one observes that the differential 
operator of 6.2 is sine: different from that of 3.1. 
The remarks which are found at the conclusion of Chapter 
III suggest the reason for this difference. The additional 
term is combined with the differential operator in order to 
indicate the improvement which follows from such an 


association. If this were not done, then 


ee iL 


and by close inspection of the results which are obtained 
below, it is evident that the sharpness of the error bound 
is materially diminished. 


Following the techniques of the earlier chapters we 


let 
a ee 
6.5 Y) = 5 


and then calculate 
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R= Ch +i r-)7t 
o 
6.6 G = “35 re 
L = 0 


Even though ( satisfies the boundary conditions on § and 
for large A approximates the solution of the differential 
equation, it is evident that ~ does not closely coincide 
with w throughout V since on D(o) the approximation is very 
poor. ; 
The set of auxiliary functions which were introduced 


in the remarks accompanying equation 3.23 may be defined as 


follows: 
Be? fi = xi i =1 and 2 


Thus we have 


6.8 ee My, =O 
M> m3 ey Mp = Q. 


In what follows the solution of 6.2 is bounded at 


the point p = (0,1). From 3.16 we have 


6.9 lw (0,1) ~2( )18| Law ne ¥T"(%)dv| 
4?! V(x) 
i, iM ¥p if (x, yw, vw) - Fy, ¥ 7) tay | 


+1 Sf vrs ttf SCP adsh* 
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where 

6.10 u, = [a (1-y)] “> : exp |- -r? ((1-y)) ~}] 

and io 

6.11 ¥/ y = ae vw =F _& 
p 4 ay 


Since Lp is such that & is zero on S, we find it 
advantageous to choose the arbitrary constants of 3.35 in 
a slightly different manner than in equations 3.38 to 3.2. 


Thus by setting 


6.12 a, rae: 

and 

6. Q =0, 
13 h 


equation 3.35 may be rewritten as 


6.1) pene jAv + (Kt- eB) SI Edy 
SZ Lp (2) Jha Aas 


2 J (A + £0) "ay aK ( I--4)| dv 


Substitution of this into 3.2) yields 





L3G 
rm Cr 
6.15 IY See Op (2) t roe {1+ expo (2K) t 


t-din -g, + BE 4 As 7 *) SSP "dv 


Due to the form of the above we set 
Xx / A 

G.10 kK=-/ # Fat +- i Ff 
a ne 


and then seek to determine the positive constants Ls and 
A. so that the right side of 6.15 is minimized. To solve 
such a@ minimum problem is extremely difficult; hence in 
practice we ordinarily choose the arbitrary constants in 
Some admissible manner and then seek to improve our result 


by a better choice of the approximation function . For 


example we assume that 


Oa? oo Sag 
16 
and set 
6.18 Soe ahi 
and 
= L we eee 
6.19 A,=R(I- Fe /6 £) a B Ar 


From the above and 6.16 it follows that 


6.20 K = 0. 





1h0 


Equation 6.15 may now be rewritten as 


sf 
6.21 ves < q+ ARI (F~E) = G. 


We use 6.10 to obtain 
Ls Le 
6.22 1S (Arsh t= {.157* = 384 


and hence 


A {SS u*ds T* {SSC By ds }* 


cy 
¢ ay R* 


where 8 is given by 6.21. 

To bound the first term on the right of 6.9 we 
introduce an auxiliary function Gr anaicn is such that the 
conclusions of Theorem 2.2 are valid. The auxiliary 
function used in this example is not identical to that 
defined by 2.39. While the original function could be 
used for the present problem, advantage is made of the 
simplicity of the differential operator to obtain a function 
which lends itself to the integrations which are required 
for the determination of the desired bound. Let 


6.2h ore (ley) exp -r@(h (1- -y)y"| 


Clearly 





11 
6.25 Th 20 for (x,y)€éV 


and since 


‘F UF | OF 
-_ TOF) = Fe + Set Sy 


=~ (ng) Benp bP Gordy] 


for (x,y) EV, it follows that the conclusions of Theorem 
2.2 are valid as was indicated above. For the present 


differential operator equation 3.57 becomes 


6.27 Wee iM {2% VI(¢) - DT@) 
Nl V4) 
a 


an x 
F>e ra(nenF Oddy 


> 


=-S(h 3F + {FO op) als. 


From 302 and the above we have 
L 
Se > Paty awa 
6.28 K K>5 > ( L + 2) <0 


where the strict inequality holds on the right due to 6.17. 


In what follows we set 


6.29 Ky = 0 





1h2 


to simplify the numerical calculations. (If greater 


accuracy were desired, one would of course take Ky equal 


to a -l}. From equations 3.57 to 3.65 we obtain 
6.30 SAF YE-v* F&)} dv 


< Bo mon {- 5} ee) Ye 


To use 3.69 it is also necessary to compute 

ena J} {(-3@ FF ep ak, (1-g)] ol V 
=e SJ (I--4) 4g 4 £2) byp lan, Ci-g)-Aec-g)y | dV 
<2. SS Ug) A tengo lntla Cin) | dv 


= 324A" 


and 


6.32 STI Pedy 


i} 


Cay S12 Cima) ® aug Fat (9) dv 


On a7 


By combining the foregoing and equation 3.69, it follows 


that 





6.33 


14 3 


|Lims SS ¥ T*(xp) dv| 
a V(4) 


BL 


a at -2)4 
< 0,204 (/,8F " {+207 ro|F Ee AGAR GE) Jt 3.29 } ) 


To bound the second term on the right of 6.9, we 


could use the techniques associated with equations 3.69 


to 3.753 however, since 


6. 3h 


oad Sp) ag Tos 


we apply Schwarz's inequality to the expression under con- 


sideration and write 


6.35 


Doon sp (FG, 4, wow) - #4, ¥, vu) fd | 
{i 
46 


<a SS %pl¥ldv 
cat SE Te iy "cat ere dy} * 
< orev at Crea stil S + Momemey 


: Y,. 
pene A ee a7 > ) 





Lyd 


The bound for the solution is completed by calculating 
the function 7 (p) which is given by 3.17. One easily 


obtains 
6.36 Up (0, /) + SJ L 8p ds =O 
6.37 NG Ve ds =—=oG¢ 
Dio) Pr 
and 
6. 38 —azizn® <— Sf (Fdv <-2212(Atk) * 


From 6.9 and the above we write as the desired bound for 


the solution 


6.39 —~o.3s¥{ ot mea EY I ~0,20H (I ra Xie” 
+ {4217 + 1.64[% + a= 
AeATNZ-2,\ + 3.24% A 


+0,06 +2212 (Att) ~ 


£< w(o,/ 
<40,384{ 3 + REAPER) t+ 026m (4A)CI.8 A" 
+ {4,217 + he4[F + AE-Z) +3,24A%} *) 


10.06 +. 22In ene 
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For example let us suppose that 
6.10 A=h 
then one obtains from 6.2 the following result 
6.1 -0.95 <w(0,1) < 41.09 


When bounds for the solution of more complicated 
problems are desired, the techniques are exactly the same 
as used above; however the integrals which result may be 


quite difficult to evaluate. 





10. 


is 
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